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Preface 


This  tbHlt  describes  as  1  avast  lgat  ion  la  to  ons  aspect 
of  s  bow  application  of  lasers  in  tbo  field  of  ocular  Medi¬ 
cine— tbe  uso  of  a  NdiYAG  laser  to  reaove  opacified  posterior 
capsules  fros  lens- lap last  patients.  Vltb  this  tool,  tbe 
opbtbalaologlst  bas  tbe  capability  of  treating  thousands  of 
people  In  tbe  office  rather  than  in  tbe  operating  room.  Tbe 
benefits  of  this  type  of  procedure,  and  a  keen  interest  in 
tbe  field  of  aediclne  were  tbe  prlae  soil rating  factors  lo 
ay  dsclsion  to  explore  this  area. 

Ibis  thesis  could  not  base  bsen  completed  without  tbe 
help  and  support  of  aany  psople,  both  physically  and  aentally. 

Without  a  doubt,  a  special  note  of  thanks  oust  be  exten¬ 
ded  to  Dr.  Leno  8.  Pedrottl ,  ay  thesis  advisor,  for  bis 
guidance,  extreae  personal  intersst ,  and  overall  supporting 
attitude.  Bis  input  was  aade  all  tbe  aore  special  as  l  was 
tbe  last  in  a  long  line  of  bis  thesis  students.  Also,  I 
wish  to  thank  Lt .  Col.  Villiaa  r.  Bailey,  who  assuasd  tbe 
duties  of  thesis  advisor  upon  the  retireaeot  of  Dr.  Pedrottl, 
and  who  also  steered  ay  efforts  in  tbe  right  direction  on 
aore  than  one  occasion. 

I  also  wish  to  thank  Majj.  Michael  B.  Staaa  for  hie 
thoughtful  questions,  and  a  grant  dsel  of  assistance.  With¬ 
out  bis,  the  STS880C  code  could  never  have  been  written. 

A  very  laportant  pert  of  tbe  work  done  in  this  project 
wen  tbe  actual  aeasureaent  of  tbe  overpressures  in  a  test 


s 
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e«U.  Tbit  could  never  have  been  accomplished  without  the 
technical  expertise  of  Boa  Gabriel  and  George  Gergal ,  both 
of  whoa  I  aa  deeply  Indebted  to. 

Thanks  aunt  also  go  to  Or.  Richard  Keates,  Jill  McCune, 
and  the  rest  of  the  etaff  at  the  Ohio  State  University  Bye 
Cllnlci  without  their  hospitality  and  help,  auch  of  this 
worh  could  never  have  been  atteapted. 

Finally,  1  wish  to  thank  Diane  for  being  so  supportive. 


keeping  as  going  when  it  » 


lapossible,  and  for  just 


being  there  when  1  truly  needed  her  the  aost. 


John  F.  Turner 


(This  thesis  was  typed  by  Its.  Barbara  A.  Grahaa. ) 
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Abstract 

Tbe  use  of  s  neodymium :YAG  laser  to  open  opacified 
posterior  capsules  in  the  human  eye  was  examined.  Hha — 
particular ,  this  study  involved  the  development  of  a  one- 
dimensional  hydrodynamic  code  to  model  the  propagation  of 
tbe  shock  wave  in  the  vitreous ,  and  the  measurement  of 
overpressures  induced  in  a  model  system. 

Tbe  model  system  consisted  of  balanced  salt  solution 
in  a  cubical  test  cell.  A  pressure  transducer  was  affixed 
to  the  rear  of  the  cell,  and  15  nsec  Q-swltched  laser  pulses 
were  focused  at  distances  of  6,  8,  10,  and  14  mm  from  the 
transducer  to  initiate  shock  waves.  Overpressures  were 
measured  for  pulse  energies  of  3.75,  7.50,  and  22.50  mJ  at 
tbe  focal  point.  V  The  balanced  salt  solution  was  chosen 
because  It  resettbles  tbe  vitreous  in  breakdown  characteris¬ 
tics,  and  the  14  mo  distance  Is  significant  since  it 
approximates  the  distance  from  tbe  posterior  capsule  to  the 
retina. 

»  Overpressures  ranging  from  0.018  -  0.432  atm  were 
recorded  at  the  transducer.  The  code  was  tested  for  a  pulse 
energy  of  3.75  mJ  and  predicted  a  pressure  change  of  2%,  in 
excellent  agreement  with  measured  values.  For  pulse  energies 
of  7.50  and  22.50  mJ,  the  maximum  deviation  was  SX  from 
measured  pressure  differences.  X 
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Based  on  current  research  by  ophthalmologists  using 
such  laser  systems  for  treatment,  it  is  concluded  that  over 
pressures  of  0.10  atm  or  less  cause  no  retinal  damage. 
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I .  Introduction 


Recent  advances  in  the  integration  of  lasers  into  the 
field  of  ophthu lmology  'mve  aided  in  the  facilitation  of 
patient  care.  One  of  these,  the  use  of  a  neodymium  YAG 
( Nd : YAG )  laser  to  remove  opacities  in  the  posterior  capsule, 
has  helped  to  move  many  patients  from  the  operating  room  to 
outpatient  cart>  facilities.  This  thesis  examines,  in  detail, 
one  consequence  of  using  a  NdsYAG  laser  to  perform  such  an 
operat ion. 

Object  I ve  and  Scope 

This  paper  will  analyse  the  shock  wave  generated  in 
the  eye  when  high  peak  power,  nanosecond  laser  pulses  are 
focused  in  the  posterior  chamber  of  the  human  eye.  The 
formation  and  propagation  of  shock  waves  in  the  posterior 
chamber  of  the  eye  ar«  modeled  and  effects  on  ocular  tissue 
are  evaluated.  The  theoretical  portion  of  this  thesis  will 
develop  a  one-dimensional  hydrodynamic  code  in  spherical 
coordinates  using  a  t.ngrangian  mesh  to  model  the  propagation 
of  a  shock  wave  through  the  eye  from  the  posterior  capsule 
to  the  retina.  The  experimental  portion  of  the  thesis  will 
involve  attaching  a  pressure  transducer  to  a  test  cell 
filled  with  balanced  salt  solution  to  measure  the  pressure 
variations  induced  by  the  pulses  of  a  Q-switched  Nd:YAG 
laser. 

The  purpose  of  this  Investigation  is  two-fold.  First, 
because  the  surgical  procedure  is  currently  being  performed 


with  Nd:YAG  lasers,  an  analysis  of  the  possible  damage  levels 
at  the  retina  due  to  the  presence  of  the  shock  wave  Is  Manda¬ 
tory.  Second,  a  correlation  between  pulse  energies  focused 


in  the  posterior  chunkier  and  actual  pressure  increases  at  the 
retina  is  highly  desirable.  This  thesis  will  develop  a 
Maximum  level  or  acceptable  overpressure  at  the  retina  during 
Nd:YAG  laser-assisted  posterior  capsulotomies.  This  maximum 
level  will  be  based  on  predicted  results,  measured  values, 
and  current  research  by  practicing  ophthalmologists. 

Background 

The  procedure  being  investigated  has  developed  from 
studies  on  the  alleviation  of  post  cataract  surgery  compli¬ 
cations.  In  such  cataract  operations,  a  surgeon  often 
removes  the  occluded  lens  from  the  patient's  eye,  and  the 
lens  is  replaced  by  external  spectacles  or  contact  lenses. 
However,  increasing  numbers  of  ophthalmologists  are  now 
UMing  intraocular  lens  (IOL)  Implants,  after  cataract 
removal ,  in  place  of  or  in  addition  to,  spectacles  or 
contact  lenses.  There  are  currently  two  methods  of  lens 
implantation  in  use,  the  posterior  chamber  method,  and  the 
anterior  chamber  method.  The  posterior  chamber  method,  the 
most  common  of  the  two  methods,  involves  placing  the  IOL 
in  the  area  vacated  by  the  patient's  own  lens  (Ref  1:271). 

The  anterior  chamber  method  of  Implantation  places  the  IOL 
anterior  to  the  original  lens  cavity,  and  requires  either 
clipping  the  1 01.  to  the  iris  (iris  clip  IOL),  or  allowing 


extensions  into  the  posterior  chamber  to  hold  the  implant 
In  position  ( iridocapsular  IOL)  (Ref  2:58). 

For  many  lens  implant  recipients,  a  serious  complication 
may  ariso.  The  posterior  capsule.  Immediately  behind  the 
patient's  lens,  opacifies,  and  reduces  the  patient's  vision 
greatly,  for  both  types  of  10L  implants.  Until  recently, 
only  further  surgery  could  remove  the  opacified  capsule  and 
restore  some  degree  of  normal  vision  to  the  patient.  However, 
in  the  past  two  years,  work  began  on  the  use  of  a  laser  to 
open  the  capsule.  By  using  a  laser,  no  surgery  in  which  the 
patient  requires  hospitalisation  is  necessary,  the  extent 
of  patient  discomfort  is  minimal,  and  the  laser-assisted 
procedure  is  accomplished  virtually  within  minutes. 

In  1980,  Aron-Rosa  of  France  began  a  study  using  a 
laser  system  to  open  the  posterior  capsule  of  lens  Implant 
patients  (Ref  3).  This  system  Is  a  mode-locked  Nd:YAG  laser 
which  emits  30  picosecond  pulses  enclosed  in  an  envelope 
whoso  pulsewldth  is  20  nanoseconds,  and  is  operated  at  a 
one  Hertz  rep  rale.  The  laser  is  focused  to  a  spot  of  50 
microns  in  diameter,  and  yields  energies  per  pulse  of  3  to 
4.5  milliJoules.  These  values  convert  to  power  densities 
of  1  x  10* 7  to  1.0  x  10* 7  I/cm7  at  the  focus. 

The  results  obtained  by  Aron-Rosa  were  very  significant. 
Twenty  post  lens-implant  patients  were  subjected  to  the  laser 
procedure  for  opening  the  posterior  capsule.  A  series  of 
five  to  ten  laser  pulses  were  focused  on  the  capsule,  and 
no  deleterious  effects  were  noted,  including  no  I0L  implants 
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being  shattered.  Effects  such  as  small  Markings  and 
scratches  were  seen  on  the  lMplanta,  however  these  would 
in  no  way  interfere  with  a  patient's  vision.  Specifically, 
no  permanent  intraocular  pressure  rise  was  noted  due  to  the 
laser  Induced  shock  wave,  and  no  lesions  were  induced  on 
either  the  cornea,  iris,  vitreous,  or  retina.  Most 
importantly,  visual  Improvement  was  achieved  in  each  patient, 
with  visual  improvement  as  great  as  4/200  to  20/20  being 
noted. 

In  1981,  Fankhauser,  van  der  Zypen,  and  associates 
used  a  similar  system  in  Germany  to  investigate  ocular 
treatm€>ntM  closely  related  to  the  opening  of  the  posterior 
capsule.  The  operations  described  by  Fankhuaser,  et.  al . . 
included  opening  pupillary  membranes  and  operations  on  the 
iris  (Ref  4). 

In  early  1982,  Lt.  James  Riggins  of  the  Air  Force 
Institute  or  Technology,  conducted  a  feasibility  study  under 
the  direction  of  Dr.  Richard  Keates  of  the  Ohio  State 
University  Eye  Clinic  and  Dr.  Leno  S.  Pedrottl  of  the 
Air  Force  Institute  of  Technology,  investigating  the  use  of 
a  g-Hwitched  laser  system  to  open  opacified  posterior 
capsules  (Kef  5).  In  addition,  the  study  investigated  the 
use  of  the  same  laser  system  to  break  vitreous  strands,  or 
cellular  striations  present  in  the  vitreous  humor. 

This  system  delivered  a  maximum  pulse  energy  of  10.89 
milliJoules  in  a  pulsewidth  of  27  nanoseconds  focused  in  a 
spot  of  diameter  4.4  microns  to  yield  a  peak  focused  power 
density  of  2.65  x  10I?  W/cm7.  Experimental  results  indicated 


that  the  system  could  open  opacified  capsules,  and  also 
that  glass  IOL  implants  would  be  shattered  it  the  IOL  were 
In  direct  contact  with  the  capsule  during  the  procedure. 
Theoretical  calculations  also  led  to  the  conclusion  that 
vitreous  strands  located  no  further  than  4-6  an  from  the 
posterior  capsule  could  be  ruptured  in  an  undiluted  eye  — 
an  eye  whose  pupil  is  the  smallest,  and  therefore  represents 
the  worst  cuse--without  possible  daaage  to  the  retina  froa 
laser  energy  depost  ion. 

The  work  previously  accoapllshed  clearly  indicates 
that  a  Nd : Y AG  laser  can  be  used  to  open  an  opacified  poster¬ 
ior  capsule  in  u  lens  iaptant  patient.  The  aechanlsa  Is  the 
creation  or  a  spark  plasm  and  an  associated  shock  wave. 

It  is  the  characterisation  of  this  shock  wave  which  is 
needed  at  this  Juncture  in  the  research.  According  to 
Dr.  Sidney  Lerman,  Professor  of  Ophthalmology ,  Emory 
University, 

"I  would  strongly  recommend  such  an 
approach;  l.e.  pressure  measurements  followed 
by  spectroscopic  and  histopathologic  evaluation 
particularly  of  the  vitreous  with  respect  to 
the  rorroer  and  the  retina  with  respect  to  the 
latter  in  order  to  come  to  some  reasonable 
conclusion  regarding  the  safety  of  this  form 
of  therapy  for  general  clinical  use."  (Ref  6:6) 

This  thesis  is  directed  toward  a  better  understanding 


of  the  generation  and  propagation  of  the  shock  wave,  and  its 
effect  on  critical  parts  of  the  eye,  especially  the  retina. 


Development 

This  thesis  will  integrate  a  portion  of  the  study  of 
hydrodynamics  with  the  basic  anatomy  and  physiology  of  the 
human  eye.  Chapter  II  will  present  the  physiology  of  the  eye 
and  the  effects  of  a  laser  on  ocular  tissue.  Chapter  III 
will  deu I  with  shock  waves,  discussing  their  generation  and 
propagation,  and  possible  effects  in  the  eye. 

Chapter  IV  will  begin  the  major  thrust  of  the  thesis, 
with  a  discussion  of  the  hydrodynamic  code  written  to  model 
shock  propagation  in  the  eye.  Chapter  V  will  explain  the 
actual  experimental  systems  and  procedures  used  to  measure 
pressure  levels  aTter  shock  generation,  and  Chapter  VI  will 
present  the  measured  results  and  evaluate  their  validity 
based  upon  predicted  data.  Finally  Chapter  VII  will  give 
clinical  conclusions  and  recomnendat ions  for  additional 
research  and  follow-up  studies. 
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II .  The  Eye  and  Laser  Radiation 


Introduction 


The  human  eye  i h  a  very  coaiplex  optical  system  designed 
to  provide  visual  sensory  input  to  the  brain.  It  is  perhaps 
the  most  important  organ  dedicated  to  the  senses,  and  any 
discussion  of  the  eye,  or  processes  occurlng  therein,  must 
begin  with  a  detailed  examination  of  the  eye  itself.  This 
chapter  begins  with  a  discussion  of  the  anatomy  and  physio¬ 
logy  or  the  eye,  then  presents  the  effects  of  laser  radiation 
on  ocular  tissue,  from  both  mechanical  and  radiative  view¬ 
points.  With  a  firm  understanding  of  the  underlying  physio¬ 
logy  of  the  eye,  a  discussion  of  the  physics  affecting  the 
organ  can  then  proceed. 


Anatomy  and 


iology  of  the  Human  Rye 


id 


Prom  an  optical  standpoint,  the  human  eye  is  essen¬ 
tially  a  positive  lens  system  which  collects  light  rays  and 
focuses  them  to  form  a  real  image  on  a  light  sensitive 
surface  (Ref  7:138-141).  Figure  1  shows  the  basic  anatomi¬ 
cal  arrangement  or  the  eye.  The  eye  is  almost  spherical, 
and  is  housed  in  a  tough  shell  called  the  sclera.  The 
sclera  is  white  and  opaque  to  light  except  for  the  front 
portion,  the  cornea ,  which  is  transparent.  It  is  the  cornea 
which  serves  as  the  first  element  of  the  lens  system.  The 
curvature  of  the  cornea  serves  to  direct  light  entering  the 
eye  through  the  anterior  chamber  which  contains  a  clear 
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watery  fluid  called  the  aqueous  huaor.  Contained  within  the 
aqueous  is  the  portion  of  the  eye  —  the  iris  —  which  gives 
the  eye  its  color.  The  iris  is  a  diaphragm  which  controls 
the  aawunt  of  light  entering  the  eye  by  expansion  or  con¬ 
traction  of  the  pupil.  The  iris  can  close  the  pupil  in 
bright  light  to  approxiaately  2  an  in  diaaeter,  or  expand 
it  in  darkness  to  about  8  on  in  diaaeter.  Also,  the  iris 
contracts  to  increase  inage  sharpness  when  the  eye  is  used 
for  close  work. 
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Figure  1.  The  human  Eye  (Ref  8:32) 
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Immediately  behind  the  iris  is  the  crystsiiins  lens, 
the  lens  is  enclosed  in  sn  envelope  --  the  so -cal led  lens 
capsule  —  which  has  walls  approximately  10  microns  thick. 
The  capsule  segment  racing  the  iris  is  the  anterior  capsule, 
while  the  segment  to  the  rear  of  the  lens  Ib  the  posterior 
capsule.  The  lens  itself  has  the  size  and  shape  of  a  small 
bean,  and  is  formed  of  approximately  22,000  layers.  The 
lens  is  initially  very  pliable}  however  this  pliability 
decreases  as  the  eye  ages.  Posterior  to  the  lens  and  its 
capsule  is  another  chamber  filled  with  a  transparent  sub¬ 
stance  composed  mostly  of  water,  the  vitreous  humor,  the 
vitreous  provides  support  for  the  retina  while  allowing 
light  and  nutrients  to  diffuse  from  the  anterior  portion  of 
the  eye  to  the  posterior  portion  (Ref  9:252). 

Within  the  sclera  is  an  inner  layer  called  the  choroid. 
This  is  a  dark  shell  containing  blood  vessels  and  melanin, 
a  dark  brown  or  black  pigment.  It  is  this  pigment  which 
k1v«*h  the  choroid  its  property  of  acting  as  an  absorber  of 
stray  light  rays  in  the  eye.  Near  the  back  of  the  eye,  a 
layer  of  light  receptor  cells  covers  much  of  the  inner 
surface  of  the  choroid  to  a  depth  of  approximately  0.51  mm. 
This  layer  is  known  us  the  retina.  A  focused  beam  of  light 
is  absorbed  in  this  structure  through  electrochemical 
reactions  occuring  in  the  two  types  of  photoreceptor  cells 
present,  rods  and  cones.  Approximately  125  million  of  these 
cells  are  contained  heterogeneously  in  the  retina.  The  rod 
cells  are  highly  sensitive  and  can  respond  to  very  dim 
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light,  however  these  ceils  do  not  distinguish  colors  and 
relay  only  ill -defined  laages.  The  cone  cells,  of  which 
there  arc  6-7  million,  give  detailed  color  laages  in  bright 
light,  but  do  not  respond  to  iaages  at  lower  light  levels. 

Contained  within  the  retina  is  the  point  of  exit  of  the 
optic  nerve  froa  the  eye.  This  region  contains  no  receptor 
cells  and  therefore  does  not  contribute  to  the  visual 
process.  This  area  is  called  the  blind  spot  for  this  reason. 
The  optic  nerve  enters  the  eye  at  the  blind  spot  and  spreads 
out  over  the  bark  of  the  interior  of  the  eye.  This  exten¬ 
sion  of  the  optic  nerve  is  the  retina  itself. 

Just  as  there  is  a  portion  of  the  retina  which  offers 
no  receptor  cells,  so  there  is  a  portion  in  which  the  cones 
are  thinner  and  more  densely  packed  than  anywhere  else. 

This  is  a  tiny  rod-free  region  approxiaately  0.3  aa  in 
diaaeter  in  its  center  known  as  the  fovea  centralis.  The 
rovea  lies  at  the  center  of  the  yellow  spot  or  aacula,  a 
saw  11  depression  2.5-3  aa  in  diaaeter  at  the  center  of  the 
retina.  The  aake-up  of  the  fovea  allows  it  to  provide  the 
sharpest  and  aost  detailed  visual  inforaatlon.  This  causes 
the  eye  to  aove  continuously  so  that  the  light  coaing  froa 
the  object  to  be  seen  is  directed  onto  the  fovea.  However, 
the  iaage  directed  onto  the  fovea  is  not  kept  stationary  on 
one  set  of  photoreceptors,  but  rather  is  constantly  shifted 
across  different,  cells.  If  the  iaage  were  stationary,  it 
would  eventually  fade  out  as  the  nerve  cells  to  which  the 
photoreceptors  are  attached  reached  saturation  —  due  to 


over-vtiaulation  —  and  canned  to  roapond.  This  la  aaialy 
due  to  the  fact  that  the  cones  in  the  fovea  are  individually 
connected  to  nerve  fibera,  while  in  contraat,  rode  are 
Multiply  connected  to  nerve  fibera,  and  any  one  fiber  can  be 
atiMulated  by  any  one  of  approxlaately  one  hundred  rod  cel la 
The  ablfttng  of  iaagea  acroaa  the  fovea  reaulta  in  the  fact 
that  the  vtaual (nation  of  any  object  ia  a  conet ruction  of 
the  ocular  ay8teM  baaed  on  a  contlnuoua  analyaia  of  a  tlae- 
varylng  retinal  iaagc. 

Figure  2  abowa  the  lualnoua  efficiency  of  the  eye. 

Both  pho topic  (adapted  for  vlalon  in  bright  light),  ?,  aad 
acotopic  (adapted  for  viaion  in  dark  light),  V,  lualnoua 


Figure  2.  Photopic  and  Scotopic  Lualnoua  Efficiency 
Curve*  for  the  Huaan  Bye  (Bef  9:362) 


efficiency  curves  are  shown,  the  aost  efficient  wavelength 
in  a  light-adapted  eye  is  555  rua ,  while  in  a  dark  adapted 
eye  it  is  510  na.  The  efficiency  of  all  other  wavelengths 
is  given  as  a  decianl  fraction  of  the  aost  efficient  wave¬ 
length  (Ref  0:352).  The  noraai  wavelength  range  of  huaan 
vision  is  approxiaaiely  froa  300  na  to  780  na.  However, 
these  liaits  have  been  extended  to  about  310  na  in  the 
ultraviolet  and  1050  na  in  the  infrared  by  soa»  studies 
(Ref  5: M0).  The  ultraviolet  transaission  liait  is  deter- 
ained  by  the  lens,  which  absorbs  in  the  UV.  Patients  such 
as  those  with  whoa  this  thesis  is  concerned  —  who  have 
undergone  surgical  lens  reaoval  —  experience  greatly 
iaproved  UV  sensitivity. 

Laser  Effect*  on  Ocular  Tissue 

Two  aaln  areas  of  specific  effects  result  when  the 
energy  froa  a  laser  is  deposited  in  the  eye.  These  effects 
can  he  classified  as  radiative  and  aechanlcal  effects. 
Radiative  effects  essentially  involve  theraal  and/or  ptioto- 
cheaical  changes  following  energy  deposition  while  aechanlcal 
duvets  deal  with  daaage  resulting  froa  acoustical  waves  or 
hydrodynaaic  shock  waves.  Even  though  only  aechanlcal 
effects  are  involved  in  the  disruption  of  opacified  poster¬ 
ior  capsules,  both  types  of  effects  will  be  exaained  in 
this  chapter. 


Radiative  Effects 


Ah  stated,  radiative  effects  involve  both  thermal  and 
photochemical  damage  mechanisms.  The  most  thoroughly 
studied,  and  adequately  modeled,  of  these  two  damage  mech - 
aniMBM  is  that  associated  with  thermal  effects.  This  type 
of  injury  accompanies  energy  absorption  in  the  retina  and 
choroid  after  foruMtng  by  the  lens  system  of  the  eye  (Ref 
MtHi).  The  heat  generated  as  a  result  of  energy  absorption 
has  two  means  or  leaving  the  focal  area  -  diffusion  and. 
conduction  by  blood  flow.  If  the  rate  of  heat  generation 
exceed**  the  rate  of  heat  conduction,  the  temperature  of  the 
tissue  rises.  If  the  temperature  exceeds  the  biologic 
tolerances  of  the  tissue,  damage  manifests  itself  as 
burning  and/or  lesion  formation. 

This  type  of  control  on  thermal  damage  has  been 
studied  by  van  der  Zypen ,  Fankhauser,  and  associates  (Ref 
1 1 : 59-180).  They  argue  against  an  increase  in  temperature 
above  the  boiling  point  in  the  iris  —  under  conditions 
of  continuous  exposure  - -  due  to  the  fact  that  local 
evaporation  initiates  fluid  circulation  in  both  anterior 
and  posterior  chambers,  whereby  hot  fluids  leave  the 
absorption  region  while  cold  fluids  enter.  This  gives  rise 
to  a  mechanism  which  limits  the  temperature  of  the  absorbing 
area.  If  a  vapor  layer  forms  between  the  absorbing  parti¬ 
cles  and  the  conducting  fluids  however,  the  temperature 
limiting  mechanism  falls  and  temperatures  above  the  boiling 
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point  may  be  reached.  Theoretical  predictions  based  on  a 
water  model  and  subsequent  experiments  conducted  using  water 
as  the  circulating  fluid  have  verified  the  accuracy  of  the 
temperature  limiting  mechanism  (Ref  11:88).  The  experiments 
were  conducted  using  a  continuous  wave  (cw)  argon  laser 
yielding  0.75  Watts  in  a  focal  spot  of  30-100  microns  in 
diamater  (Ref  11:89-90). 

Either  general  damage  mechanism  —  radiative  or 
mechanical  —  is  highly  dependent  upon  the  time  of  exposure. 
This  is  due  to  the  fact  that  for  longer  exposures,  more  of 
the  energy  is  absorbed  and  interacts  with  the  tissues  under¬ 
going  irradiation.  If  exposures  of  more  than  a  few  seconds 
are  being  considered,  then  evidence  indicates  that  a  damage 
mechanism  referred  to  as  "photochemical"  becomes  dominant 
(Rer  10:84).  These  effects  are  primarily  associated  with 
blue  and  ultraviolet  light,  and  involve  light  whose  photon 
energies  correspond  to  the  binding  energies  of  cell  consti¬ 
tuents  resulting  in  the  breaking  of  molecular  bonds  or  the 
initiation  of  chemical  reactions  in  the  tissue  surrounding 
the  focal  spot  (Ref  12:1269). 

Photochemical  damage  manifests  itself  in  various  ways 
depending  on  the  ocular  tissue  involved.  The  damage  mecha¬ 
nism  in  the  retina  is  the  formation  of  a  lesion,  while 
damage  in  the  cornea  involves  destruction  of  epithelial 
cells,  or  in  some  more  severe  cases,  corneal  opacification 
(Ref  13:108-110).  Photochemical  damage  to  the  lens,  which 
may  require  a  great  deal  of  time  to  appear,  depending  on 
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the  dose  involved,  manifests  itself  as  the  destruction  and 
degradation  of  lens  materials,  leading  to  permanent  opaci¬ 
fications  of  the  lens  (Ref  13:106,  112-114).  Such  lens 
opacifications  are  commonly  referred  to  as  cataracts. 

Mechanical  Ef fects 

As  exposure  times  become  shorter,  mechanical  effects 
may  become  significant  and  lead  to  damage.  However,  much 
controversy  still  surrounds  the  extent  of  this  type  of 
damage.  Two  sources  of  mechanical  damage  exist,  acoustical 
waves  and  shock  waves.  Both  are  caused  by  the  absorption  of 
a  great  deal  of  energy  in  a  very  short  time  period.  However 
the  shock  wave  propagates  at  a  velocity  greater  than  the 
local  speed  of  sound  while  the  acoustical  wave  propagates 
at  a  velocity  equal  to  the  local  sound  speed  (Ref  14:111). 

Extensive  research  by  Cleary  and  associates  (Ref  15:175- 
219)  on  acoustical  wave  propagation  and  damage  has  Indicated 
that  damage  due  to  these  transients  can  be  expected  over  a 
much  larger  area  than  would  be  noticed  from  thermal  damage 
alone.  Also,  it  must  be  noted  that  while  the  two  types  of 
pressure  transients  under  discussion  are  normally  considered 
as  separate  quantities,  a  shock  wave  will  degenerate  into 
an  acoustical  wave  if  the  density  of  the  surrounding  medium 
is  sufficiently  high,  or  if  the  time  interval  of  propagation 
is  sufficiently  long.  Research  conducted  by  van  der  Zypen , 
Fankhauser,  and  associates  has  indicated  that  shock  wave 
formation  in  the  iris  can  create  hemorrhaging,  cell  f ragmen- 


tation,  and  the  ejection  of  blood  cells  into  the  surrounding 
tissue  (Ref  11:135-136).  A  similar  type  of  damage  is  seen 
on  the  retina. 

The  above  discussion  of  the  anatomy  of  the  eye  and  its 
reactions  to  laser  radiation  has  laid  the  necessary  ground¬ 
work  for  this  study.  Now  we  must  examine  the  physics  of 
shock  waves  before  we  can  fully  understand  the  processes 
which  occur  when  a  shock  wave  is  created  in  the  vitreous 
humor  or  the  posterior  capsule.  Chapter  III  provides  this 


information 


III.  The  Physics  of  Shock  Waves 

Introduct ion 

At  this  point,  we  have  gained  an  understanding  of  the 
human  eye  as  an  organ  dedicated  to  providing  visual  sensory 
input  to  the  brain.  However,  when  the  energy  from  a  Nd:YAG 
laser  pulse  is  deposited  in  the  posterior  capsule  or  adja¬ 
cent  vitreal  tissue,  the  main  processes  of  interest  become 
more  physical,  and  less  biological,  in  nature.  This  chapter 
will  provide  the  material  necessary  to  understand  the  phy¬ 
sical  nature  of  a  shock  wave,  both  in  terms  of  its  generation 
—  how  a  shock  wave  is  created  —  and  its  propagation  —  how 
a  shock  wave  moves  through  the  surrounding  medium.  In 
addition,  the  possible  effects  caused  by  a  shock  wave  in 
various  parts  of  the  eye  will  be  discussed. 

Shock  Waves 

In  the  following  discussion,  the  terms  fluid  and  gas 
will  often  be  used  interchangeably  to  denote  the  medium  in 
which  the  shock  is  propagating.  This  is  often  encountered 
in  fluid  dynamics,  where  a  fluid  is  considered  any  substance 
which  flows,  and  a  gas  is  a  subset  of  the  overall  class  of 
fluids. 

A  discussion  of  shock  waves  is  best  begun  by  examining 
a  gas  initially  at  rest,  with  constant  density  pQ  and 
pressure  pQ ,  bounded  on  the  left  by  a  plane  piston.  Assume 
now  that  the  piston  moves  to  the  right  with  velocity  u, 


such  that  the  gas  is  compressed.  Three  basic  gas  dynamic 
equations  govern  the  motion  of  the  gas  either  prior  to  the 
motion  of  the  piston  or  after  compression  begins.  These  are 
the  equations  of  continuity,  the  equation  of  momentum  con¬ 
servation,  and  the  equation  of  energy  conservation  of  a 
small  group  of  fluid  particles  commonly  located  near  a 
point  x  at  a  time  t  (Ref  16:104). 

dp/dt  *  -p(V-u) 

du/dt  *  -p-1(Vp) 

dE/dt  =  -pp-1(7*u) 

respectively,  where, 

u  =  average  velocity  (cm/s ec) 
p  =  average  mass  density  (gm/cm3) 
p  =  average  pressure  (erg/cm3) 

E  =  average  specific  internal  energy  (erg/gm) 

If  an  attempt  is  made  to  find  a  continuous  solution  for  the 
case  of  a  piston  with  no  initial  velocity  such  that  the 
velocity  of  the  piston  is  discontinuous,  no  physically  real¬ 
istic  solution  is  available.  However,  a  solution  can  be 
found  if  it  is  assumed  that  the  variables  u,  p,  and  p  are 
aiscont inuous  across  a  boundary  which  defines  regions  where 
the  gas  is  undisturbed  —  u  *  0,  p  *  pQ,  p  =  pQ  —  and  where 
the  gas  immediately  adjacent  to  the  piston  moves  with  velo¬ 
city  equal  to  the  piston  velocity  (Ref  17:45).  Such  a 
solution  is  shown  in  Figure  3. 


(1) 

(2) 

(3) 


In  this  type  of  solution  to  the  gas  dynamic  equations, 
the  transition  from  undisturbed  gas  to  compressed  gas  occurs 
almost  instantaneously  over  such  a  thin  region  that  the  dis¬ 
continuous  model  is  in  many  cases  quite  valid.  This  type 
of  discontinuity  is  known  as  a  weak  shock.  While  Equations 
(1)  -  (3)  for  an  ideal  fluid  predict  the  development  of  a 
shock,  they  become  meaningless  at  a  shock  and  cannot  be 
used  for  its  complete  treatment.  If  the  viscosity  of  the 
medium  is  introduced  in  the  equations,  the  overall  discon¬ 
tinuity  of  the  solution  relaxes  to  a  point  where  the  equa¬ 
tions  can  be  solved  to  predict  shock-transition  details 
(Kef  18:35). 


Figure  3.  Initial  State,  Density,  and  Velocity 
Profiles  of  a  Shock  Wave  (Ref  17:45) 

Shock  Wave  Generation 

Various  phenomena  are  associated  with  the  generation 
of  a  shock  wave.  The  movement  of  a  piston  through  an  un¬ 
disturbed  gas  has  been  mentioned  —  and  will  be  examined  in 


detail  subsequently  —  but  this  is  not  the  only  means  of 
creating  a  shock  wave.  Another  method  of  shock  generation 
involves  a  strong  explosion,  in  which  a  great  deal  of  energy 
is  deposited  in  a  small  volume  in  a  very  short  time  interval 
(Ref  15:93-101).  It  is  this  type  of  shock,  referred  to  as 
a  detonation  or  reactive  shock,  that  is  generated  when  an 
energy  pulse  from  a  Nd:YAG  laser  is  deposited  quickly  in  the 
posterior  capsule  of  the  eye.  Several  theories  have  been 
offered  concerning  the  shock  wave  accompanying  plasma  forma¬ 
tion.  Theroized  wave  types  include  a  detonation  wave  driven 
by  plasma  absorption  and  a  radiation  transport  wave  in  which 
the  radiated  thermal  energy  from  the  plasma  acts  to  ionize 
material  behind  the  plasma  which  then  absorbs  laser  radia¬ 
tion  (Ref  19:68,  70).  Normally  such  an  absorption  wave, 
also  known  as  a  laser  supported  detonation  wave,  exists 
when  the  laser  energy  is  being  provided  by  a  continuous  beam 
rather  than  a  pulsed  beam.  This  is  due  to  the  fact  that 
with  continuous  wave  (cw)  radiation,  the  plasma  is  constant¬ 
ly  being  fed  with  energy  that  must  be  absorbed,  while  with  a 
pulsed  beam,  the  energy  input  is  not  constant,  and  the  plasma 
can  dissipate  if  the  time  between  laser  pulses  is  great 
enough . 

The  shock  waves  resulting  from  laser  induced  plasma 
formation  have  been  studied  by  several  research  scientists 
(Ref  20;  21).  Velocities  ranging  from  2  x  105  cm/sec  to 
1.6  x  106  cm/sec  and  pressures  as  high  as  250  kbar  (1  bar  * 
0.9869  atm  =  9.97  x  105  erg/cro3)  have  been  reported  using 


.■o.o. 


both  Nd:YAG  and  ruby  lasers  whose  powers  were  usually  less 
than  5MW,  but  occasionally  as- high  as  30  MW  (Ref  20:282; 
21:46).  All  of  these  values  were  obtained  using  water  as 
the  medium  in  which  the  plasma  was  formed  and  the  associated 
shock  wave  was  developed. 

Shock  Wave  Propagation 

In  order  to  gain  an  understanding  of  how  a  shock  wave 
propagates  in  a  medium,  the  example  of  the  piston-driven 
shock  will  be  examined  in  detail.  The  piston  shock  rather 
than  the  reactive  shock  will  be  examined  because  the  piston 
shock  can  be  described  by  Cartesian  coordinates  while  the 
reactive  shock  is  often  described  by  spherical  coordinates. 
The  difference  arises  in  the  fact  that  for  spherical  shock 
waves,  the  pressure  front  moves  outward  radially  from  the 
plasma,  assuming  a  1/r2  dependence  in  a  homogeneous  medium 
where  r  is  the  distance  from  the  plasma  center  to  the  leading 
edge  of  the  shock  front.  This  l/rz  dependence  is  not  present 
in  the  piston  driven  shock,  thereby  making  it  a  much  more 
easily  understood  example  for  illustrative  purposes. 

Returning  to  the  case  of  the  moving  piston,  recall  that 
the  piston  is  moving  with  a  velocity  u,  and  is  forming  a 
layer  of  gas  ahead  of  it,  as  illustrated  in  Figure  4. 

In  this  example,  following  the  treatment  of  Stamm, 

(Kef  16:114-116),  after  a  given  time  t,  a  mass  of  gas  per 
unit  area  expressed  as  PQDt  (gm/cm2)  where  D  is  the  velocity 
of  the  gas,  is  piled  up  in  front  of  the  piston.  D  must 
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Figure  4.  The  Piston  Driven  Shock  Front 

necessarily  be  greater  than  u  to  account  for  the  accumula¬ 
tion  of  gas.  The  gas  is  compressed  into  a  length  given  by 
(D-u)t  with  a  density  of  p^.  This  forms  the  compression 
region  mentioned  previously.  The  conservation  of  mass 
requires  that  the  total  uncompressed  mass  equal  the  total 
compressed  mass,  or, 

pjCD-uJt  =  pQDt  (4) 

which  is  another  form  of  the  equation  of  continuity.  In  a 
similar  manner,  we  can  obtain  the  equation  of  momentum 
conservation  by  recognizing  that  the  mass  will  acquire  a 
momentum  per  unit  area  given  by  pQDtu  (gm/sec-cm)  which  must 
equal  the  impulse  per  unit  area,  given  by  (p^  -  pQ)t,  so 
that , 


Finally,  if  the  increase  in  total  energy  of  the  compressed 
gas  is  considered,  this  must  equal  the  energy  added  to  the 
gas  by  the  action  of  the  piston,  as  dictated  by  the  equation 
of  energy  conservation.  This  means  that, 

P  jUt  -  poDt(E1-EQ  +  iu2)  (6) 

where , 

pQDt  is  the  mass  per  unit  area 

Ej-Eo  is  the  internal  energy  per  unit  mass  (erg/gm) 
£uz  is  the  kinetic  energy  per  unit  mass  (erg/gm) 

In  order  to  cast  the  three  equations  we  have  presently  into 
a  set  that  more  closely  resemble  Equations  (1)  -  (3),  we  must 
transform  the  coordinate  system  of  interest  from  the  one 
illustrated  in  Figure  4  to  one  moving  with  the  shock  front 
at  a  velocity  D.  This  is  done  by  defining  new  velocity 
variables  uq  and  Uj  such  that, 


-D 

(7) 

-(D-u) 

(8) 

whore , 

uo  =  the  velocity  of  the  gas  flowing  into  the 
discontinuity 

Uj  *  the  velocity  of  the  gas  flowing  out  of  the 
discontinui ty 

This  new  coordinate  system  is  illustrated  in  Figure  5. 


Figure  5.  The  Shock  Discontinuity 


In  this  new  coordinate  system,  the  gas  dynamic  equa¬ 
tions  —  the  equation  of  continuity,  the  equation  of  momentum 
conservation,  and  the  equation  of  energy  conservation, 
respectively  —  can  be  expressed  as, 


plul 
P1  +  plul2 
E1  +  pl/pl  +  *ul2 


P  u 
Fo  o 

P  +  P  u  2 
o  Fo  o 

Eo  +  po/po  +  iuo5 


(9) 

(10) 

(11) 


For  any  given  set  of  initial  conditions,  another  equation  is 
required  to  completely  describe  the  motion  of  the  gas.  This 
is  required  since  we  have  three  equations  and  four  unknown 
quantities  —  Pj  ,  u^,  p1 ,  and  E^.  This  is  the  equation  of 
state,  which  gives  a  relation  between  the  internal  energy  E 
and  the  pressure  p,  both  for  the  ambient  gas  (Eq,  PQ)  and 
for  the  compressed  gas  (E^,  Pj)»  For  an  ideal  gas,  the 
equation  of  state  is  given  by, 


( Y- 1 )E  ■  p/p  =  kT/m 


(12) 


where , 


k  -  Boltzmann* 8  constant  (1.38  x  10~l°  erg/°K) 

T  •  the  absolute  temperature  (°K) 
m  -  the  molecular  weight  of  the  gas  ( gm ) 

and, 

Y  -  Cp/Cy  (13) 

where  f 

Cp  ■  the  specific  heat  of  the  gas  at  constant 
pressure  (erg/gm-°K) 

Cy  ■  the  specific  heat  of  the  gas  at  constant 
volume  (  erg/gm-°K) 

At  this  point,  we  could  analyze  any  gas  or  liquid  in  which 
a  shock  wave  is  propagating,  if  we  bad  the  true  equation  of 
state,  which  often  differs  froa  '.he  ideal  gas  condition  of 
Equation  (12). 

Shock  Wave  Kf fee is  on  Ocular  Tissue 

While  many  investigations  have  been  conducted  in  an 
attempt  to  quantify  the  amount  of  damage  directly  attribu¬ 
table  to  laser  Induced  shock  waves,  very  little  definitive 
information  is  available.  This  seems  to  follow  because  in 
roost  of  the  reported  experiments,  the  system  under  study, 
e.g.  the  retina,  was  of  a  sufficiently  complex  nature  that 
it  was  not  possible  to  differentiate  between  the  effects  of 
other  sources  of  tissue  disruption  —  thermal  denaturation , 


vaporization,  and  photochemical  inactivation  —  and  effects 
due  solely  to  mechanical  dell  disruption  caused  by  laser 
induced  shock  eaves.  In  most  of  these  cases,  due  to  the 
laser  energy  densities  and  Irradiation  regions,  the  tissue 
areas  exposed  to  shock  waves  of  sufficient  amplitude  to 
cause  damage  were  in  close  proximity  to  the  site  of  laser 
absorption,  and  in  these  areas  gross  tissue  damage  was 
visible  from  other  effects  (Ref  15:201).  More  research  is 
needed  on  this  topic;  however,  in  the  area  of  damage  to  the 
retina,  some  data  has  been  obtained. 

Retinal  Effects 

There  is  no  doubt  that  the  retina  Is  susceptible  to 
damage  from  the  shock  wave  caused  by  a  laser  pulse.  Marshall 
(Her  22:97-100)  classified  this  damage  into  two  main  areas: 
tissue  fracture  induced  by  the  passage  of  a  strong  pressure 
front  —  a  shock  wave  —  and  tissue  fracture  due  to  stresses 
caused  by  local  displacement  of  the  affected  cells.  In  the 
retina,  the  first  type  of  effect  manifests  itself  as  damage 
to  the  photoreceptors  in  the  outlying  areas  of  the  retinal 
lesion,  while  the  second  type  of  effect  forms  holes  and  tears 
in  the  nerve  fibers  posterior  to  the  receptor  cells.  Marshall 
suggests  that  the  damage  to  the  outer  segments  of  the  recep¬ 
tors  is  indicative  of  a  single  shock  front  rather  than  a 
resonance  effect  in  the  tissue  layer. 

Pressure  transients  generated  during  laser  photocoagu¬ 
lation  as  well  as  pulsed  xenon  (Xe)  retinal  photocoagulation 


were  investigated  by  Fraunfelder  and  Viernstein  (Ref  23:1261- 
1270).  This  study  involved  the  implantation  of  pressure 
transducers  directly  into  various  ocular  tissues  of  both 
rabbits  and  monkeys.  A  ruby  laser  providing  700  psec  pulses 
of  0.6  J  energy  caused  overpressures  of  up  to  130  torr 
(0.171  atm)  at  the  point  of  measurement.  The  Xe  pulses 
cuused  minimal  overpressures  but  did  cause  retinal  lesions. 
The  authors  stated  that  no  evidence  of  ocular  damage  was 
detected  from  overpressures  of  10  -  40  torr  (0.013  -  0.053 
atm)  observed  as  a  result  of  shock  wave  generation  due  to 
tuner  pulse  absorption. 

Cleary  and  Hamrick  have  studied  the  problem  of  pressure 
transients  in  the  eye  extensively.  These  authors  conclude 
that  these  pressure  waves  can  definitely  damage  ocular 
tissue  at  or  near  the  site  of  absorption,  but  not  at  dis¬ 
tances  far  removed  from  the  absorption  site.  They  suggest, 
on  the  basis  or  experimental  and  theoretical  evidence,  that 
at  the  energy  densities  associated  with  threshold  retinal 
effects,  the  mu  in  cause  of  damage,  second  only  to  direct 
thermul  effects,  is  the  overpressure  caused  by  a  laser 
induced  shock  wave  (Ref  24:1-10). 

At  this  time,  the  exact  pressure  amplitude  of  a  shock 
wave  needed  to  cause  retinal  damage  has  not  been  determined. 
However,  the  extensive  laser  treatment  of  patients  by  Aron- 
Rosa  and  Funkhauser  has  not  indicated  the  existence  of  reti¬ 
nal  damage  to  human  patients,  even  after  multiple  laser 
pulse  depositions  in  the  eye,  although  both  authors  agree 


that  it  is  the  action  of  the  shock  wave  that  destroys  the 
posterior  capsule  (Ref  3;4).  This  is  not  surprising  when 
it  is  realized  that  the  plasma  initiated  by  the  laser  pulse 
has  a  radius  on  the  order  of  tens  of  microns  while  the  retina 
is  approximately  1.4  cm  from  the  focal  point.  The  small 
initiation  volume  and  the  comparitively  distant  retina 
enables  the  pressure  of  the  shock  wave  to  be  drastically 
reduced  in  amplitude  while  transversing  the  distance  from 
initiation  site  to  the  retina. 

To  see  this  more  clearly,  we  turn  to  the  treatment  of 
a  shock  wave  introduced  by  Sedov  (Ref  25:238-251).  In  this 
particular  treatment,  the  shock  is  reduced  to  a  set  of 
dimensionless  variables  which  relate  the  initial  energy  of 
the  laser  pulse,  the  ambient  pressure  of  the  medium,  the 
pressure  of  the  leading  edge  of  the  shock,  and  the  radius  of 
the  leading  edge  of  the  shock.  A  thorough  discussion  of 
these  variables  is  included  in  Appendix  F.  If  we  consider 
an  initial  laser  pulse  energy  of  10  m J ,  based  on  this  treat¬ 
ment,  the  maximum  overpressure  we  could  expect  for  a  value 
at  the  retina  would  be  0.096  atm.  This  is  based  on  the 
assumption  that  the  laser  pulse  is  focused  in  an  ideal  gas, 
with  an  ideal  gas  constant  equal  to  1.022. 

Obviously,  the  vitreous  is  not  an  ideal  gas,  however, 
we  can  use  this  value  as  a  first  approximation  to  such  over¬ 
pressures  as  we  would  expect  to  predict  and/or  measure  at 
the  retina. 

With  this  background,  we  can  now  examine  the  effects  of 


such  a  shock  wave  in  terms  of  the  overpressure  present  at 
the  retina.  The  first  step  in  the  analysis  is  the  develop¬ 
ment  of  a  hydrodynamic  code  to  predict  such  overpressures. 
Chapter  IV  presents  this  code,  and  the  theory  necessary  for 
its  derivation. 


IV.  Theoretical  Model ing  of  a  Shock  Wave 
In  the  Eye 


Introduction 

Many  engineering  problems  exist  which  deal  with  the 
effects  on  various  sytems  of  large  energy  inputs  into  small 
volumes  in  short  time  intervals.  Several  computer  codes 
exist  which  deal  with  the  resulting  shock  waves  which  form 
and  propagate  in  such  systems  (Ref  26;  27).  The  shock  wave 
which  occurs  in  the  eye  can  be  approximated  in  much  the  same 
fashion  if  a  proper  model  for  the  eye  is  used.  This  means 
that  if  the  correct  ambient  pressure,  ambient  density,  and 
geometry  for  the  eye  is  used,  and  more  importantly,  if  the 
correct  equation  of  state  is  used,  a  finite  difference  scheme 
can  be  written  which  allows  the  prediction  of  overpressure 
levels  at  the  retina  as  a  function  of  initially  deposited 
laser  pulse  energy.  In  this  chapter,  the  development  of  a 
one-dimensional  hydrodynamic  code  in  spherical  coordinates 
using  a  Lagrangian  mesh,  EYESHOC,  will  be  presented.  This 
code  is  intended  to  yield  an  estimate  of  such  overpressures. 

The  main  objective  in  the  development  of  EYESHOC  is  to 
design  a  code  which  accurately  predicts  the  overpressure 
at  the  retina  caused  by  a  laser  induced  shock  wave.  More¬ 
over,  this  must  be  accomplished  while  including  simplifying 
assumptions  which  increases  the  ease  of  using  the  code,  and 
decrease  the  required  computer  time.  A  listing  of  the  code, 
the  definitions  of  all  variables  involved,  and  a  sample  of 
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the  output  from  EYESHOC  is  contained  in  Appendix  E 


One-Dimensional  Hydrodynamics  Code 


The  following  subsections  present  the  basic  equations 
and  assumptions  used  to  develop  EYESHOC,  and  also  include  a 
discussion  of  the  Lagrangian  mesh  which  forms  the  coordinate 
system  for  the  code. 


Governing  Equations 


The  basic  equations  which  govern  the  flow  of  a  fluid 
were  presented  in  Chapter  III.  These  are  the  equation  of 
continuity,  and  the  equations  of  momentum  and  energy  conser¬ 
vation,  Equations  (1)  -  (3).  As  was  discussed,  an  additional 
relation  between  pressure  and  internal  energy,  the  equation 
of  state,  is  needed  to  complete  the  set  of  differential 
equations  required  for  the  calculation  of  particle  motions. 


Unfortunately,  the  vitreous  humor  is  not  an  ideal  gas,  so 
Equation  (12)  is  not  sufficient  for  EYESHOC  calculations. 
Instead,  the  Gruneisen  equation  of  state  is  used.  This 
equation  of  state  is  much  more  generally  applicable,  and  is 
(Ref  18:3) 


where 


P  =  PH  +  (Yg/V)(E-EH) 


PH  *  c2(vo-v)/[vo  -  S(Vo-V)] 


(14) 


(15) 


:H  =  i  [  [C(vo-v)]/[vQ  -  s(vo-v>]  ]2 

(16) 

=  y-1 
s  ' 

(17) 

V  =  1/p 

(18) 

S  =  iy 

(19) 

and 


pQ  is  the  normal  density  of  the  material  (gm/cm3) 
p  is  the  density  (gm/cm3) 

C  is  the  speed  of  sound  in  the  unshocked 
material  (cm/sec) 

Assumptions 

The  following  assumptions  were  made  in  order  to  arrive 
at  the  differential  equations  solved  by  EYESHOC: 

1.  The  local  fluid  velocity  is  dependent  upon  and 
directed  along  the  x-coordinate  in  a  Cartesian  coordinate 
system,  or  along  the  r-coordinate  in  a  spherical  or  cylin¬ 
drical  coordinate  system,  thereby  imposing  that  the  code  be 
one-dimensional  in  nature. 

2.  Energy  input  into  the  code  is  assumed  to  be  instan¬ 
taneous,  so  that  the  addition  of  energy  into  the  system  is 
handled  as  an  initial  condition.  In  addition  to  the  assump¬ 
tion  of  instantaneous  temporal  input,  the  energy  is  assumed 
to  be  deposited  at  a  geometrical  "point"  at  the  center  of  a 
spherical  coordinate  system. 
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3.  An  artificial  viscosity,  q,  can  be  introduced  into 
the  code  to  prevent  the  development  of  numerical  instabili¬ 
ties  which  develop  at  a  shock  front  due  to  the  contraction 
of  the  material .  The  introduction  of  this  variable  allows 
for  the  shock  to  spread  out,  which  allows  for  a  numerical 
solution  to  the  problem  to  exist.  This  is  accomplished  by 
replacing  p  in  Equations  (2)  and  (3)  with  p+q  where  q  is 
defined  for  a  Langrangian  mesh  as  (Ref  28:313) 

(aAr  )2 p(V*u)2  if  V*u  <  0 

q  =  (20) 

0  if  V*u  l  0 

where 

Ar  is  the  cell  size,  and 

a  is  a  constant  equal  to  approximately  one-half 
the  number  of  cells  over  which  the  shock  is 
spread 

4.  The  applicable  equation  of  state  for  the  calcula¬ 
tions  performed  by  EYESHOC  is  the  Griineisen  equation  of  state. 
This  is  the  correct  equation  because  the  vitreous  is  of 
sufficient  viscosity  to  approach  the  properties  of  a  solid, 
and  this  equation  of  state  is  applicable  to  many  solids 

(Ref  18:3).  In  order  to  obtain  a  time  derivative  of  the 
pressure  needed  to  fit  into  the  finite  difference  scheme  of 
the  calculation,  the  specific  internal  energy  is  eliminated 
from  the  equation  of  state  as  well  as  the  rest  of  the 


m 
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necessary  differential  equations.  Appendix  C  contains  the 


derivation.  Once  this  has  been  accomplished,  the  resulting 


equation  is 


dp/dt  =  - 


p  +  ( Y-l )(p+q)  -p2pQCi 


Jy(p-p„)p, 


[p-  4y(p-pq)] 


7-u  (21) 


Lagrangian  Mesh 


The  spatial  mesh  used  in  EYESHOC  is  defined  in  a 


Lagrangian  coordinate  system  such  that  each  fluid  particle 


is  labeled  with  a  Lagrangian  coordinate  r  which  is  retained 


as  the  particle  moves  about.  The  volume  between  adjacent 


mesh  points  is  defined  as  a  cell,  and  the  boundaries  of 


these  cells  are  free  to  move  about  as  necessary,  but  only  in 


the  positive  or  negative  r-directions.  These  cell  bounda¬ 


ries  form  concentric  spheres  in  the  case  of  the  EYESHOC  mesh. 


In  essence,  the  boundaries  in  a  Lagrangian  mesh  appear  as 


masses  connected  by  springs  since  they  are  free  to  move  from 


point  to  point  (Ref  16:105;  16:82;  29:14-15),  where  such 


springs  would  correspond  to  the  connecting  variables  such  as 


pressure,  density,  etc. 


The  overall  scheme  present  herein  is  illustrated  in 


Figure  6.  Figure  6a  shows  the  geometry  of  the  eye,  with  the 


point  of  laser  energy  deposition  identified,  and  the  resulting 


shock  waves  graphically  represented.  Figure  6b  shows  how 


these  shock  waves  would  propagate  in  a  Cartesian  coordinate 


system,  while  Figure  6c  gives  the  resulting  shock  wave 


propagation  converted  to  a  spherical  coordinate  system, 
thereby  delineating  the  development  of  the  Lagrangian  mesh 
used  in  EYESHOC. 


The  variables  appearing  herein  are  either  defined  at 
the  boundaries  of  the  cells  or  at  the  centers  of  the  cells. 
Particle  velocity  and  position  are  defined  at  cell  boundaries. 
Pressure,  artificial  viscosity,  density,  and  specific  internal 
energy  are  defined  at  the  cell  center.  Figure  7  shows  these 
quantities  in  the  spatial  mesh. 
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Figure  7.  Lagrangian  Cells  and  Cell  Variables 
Used  In  EYESHOC 


By  convention,  in  a  Lagrangian  system,  the  origin  ;>f  the 
radial  coordinate  is  chosen  to  be  the  position  of  the  rQ  cell 
boundary.  This  is  the  center  of  the  first  cell  at  zero  time. 


This  point  is  fixed,  and  does  not  vary  as  the  positions  of 
the  cell  boundaries  are  updated  by  each  successive  iteration 
of  the  code.  It  should  be  noted  here  that  the  cells  are  in¬ 
deed  concentric  spheres,  however  this  was  not  illustrated 
in  Figure  6  in  order  that  the  variables  and  their  relations 
to  the  cells  could  be  shown  more  easily. 

When  dealing  with  hydrodynamic  shock  codes,  an  addi¬ 
tional  coordinate  system  could  be  considered,  this  being 
the  so-called  Eulerian  mesh  (Ref  18:82).  In  this  mesh,  the 
cell  boundaries  remain  fixed  while  the  fluid  particles  move 
through  the  mesh.  The  Lagrangian  mesh  is  more  accurate 
when  dealing  with  a  free  surface,  however  it  can  become 
inaccurate  when  the  cells  become  strongly  distorted.  An 
Eulerian  mesh  allows  arbitrarily  great  distortions,  however 
an  Eulerian  mesh  is  virtually  unable  to  resolve  the  shape 
of  a  free  surface.  The  Lagrangian  mesh  was  considered  to 
be  the  most  logical  choice  for  this  thesis. 

Computer  Model  of  the  Human  Eye 

This  section  will  present  the  model  of  the  human  eye 
used  in  EYESHOC  to  predict  overpressure  levels  present  at 
the  retina.  The  limiting  assumptions  will  be  discussed  and 
the  initial  conditions  which  need  to  be  determined  for  input 
into  the  code  will  be  examined,  as  will  the  basic  scheme  of 
the  code  design. 

Assumptions 

The  following  assumptions  were  made  in  order  to  further 
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modify  a  standard  one-dimensional  hydrodynamic  code  as 
described  in  the  earlier  sections  of  this  chapter  to  model 
shock  propagation  in  the  eye: 

1.  The  eye  is  a  sphere  and  shock  waves  propagating 
within  it  will  also  be  spherical  in  nature,  varying  only 
in  the  radial  or  r-direction. 

2.  The  eye  is  filled  with  a  homogeneous  fluid  medium 
of  density  1.007  gm/cm3  (Ref  9:252)  which  has  virtually  the 
same  properties  as  water,  and  can  therefore  be  approximated 
as  such,  with  the  exception  that  the  viscosity  of  the 
vitreous  is  much  greater  than  that  of  water.  From  this 
approximation,  a  value  of  y  can  be  calculated.  Appendix  A 
provides  this  calculation,  with  the  result  that  for  water  — 
and  therefore  for  the  vitreous  also  —  y  ■  1.022. 

3.  For  the  calculations  performed  by  this  code, 
reflections  of  the  spherical  shock  waves  will  not  be  consid¬ 
ered  in  regard  to  the  overpressure  levels  at  the  retina.  At 
most,  these  reflections  will  amount  to  an  increase  in  the 
overpressure  of  a  factor  of  two  (Ref  15:189).  For  over¬ 
pressure  levels  such  as  those  reported  in  Chapter  III,  this 
is  not  considered  to  be  significant.  Reflections  from  the 
roof,  floor,  and  front  of  the  eye  will  therefore  be  ignored. 

4.  Posterior  capsulotomies  will  only  be  modeled  in 
patients  who  have  had  their  lens  removed,  and  who  have  not 
had  a  glass  or  plastic  IOL  implanted  in  their  eye.  This 
serves  to  alleviate  the  need  to  consider  shock  reflections 
from  these  substances,  as  well  as  the  effect  on  the  shock  of 


propagation  into  either  a  human  lens  or  an  IOL. 

5.  The  distance  from  the  posterior  capsule  to  the 
retina  is  taken  to  he  a  nominal  14  mm. 

6.  All  of  the  energy  from  the  laser  pulse  is  converted 
into  internal  energy  for  the  purpose  of  initiating  the  shock 
wave.  This  energy  is  deposited  instantaneously  and  the  shock 
front  develops  immediately  thereafter. 

7.  The  ambient  pressure  of  the  eye  is  780  torr,  or 
1.026  atm.  This  ambient  pressure  does  not  affect  the  nature 
und/or  propagation  of  the  shock  wave. 

Code  Discussion 

EYESHOC  was  designed  as  an  easy-to-use  analytical  tool 
to  be  employed  by  anyone  needing  the  type  of  data  it  can 
provide.  The  initial  conditions  required  of  the  user  are 
very  few,  and  will  be  discussed  subsequently.  The  important 
features  of  the  code  are  as  follows: 

1.  The  code  is  allotted  50  spherical  cells,  with  the 
origin  of  these  cells  at  the  center  of  cell  one.  Each  cell 
is  larger  than  the  previous  cell  by  a  factor  of  10%,  i.e., 
r2  equals  l.lOr^;  r3  equals  l.lOr^;  r4  equals  1.10r3,  and 
so  forth  up  to  r5Q. 

2.  Of  these  50  cells,  only  30  are  actually  needed  to 
allow  the  shock  to  travel  the  distance  from  the  posterior 
capsule  to  the  retina.  The  remaining  cells  are  included  so 
that  edge  effects  —  numerical  instabilities  which  arise 
when  the  right-most  edge  of  the  coordinate  system  is  en- 


countered  during  the  iterative  process  —  can  be  ignored.  In 
its  current  form,  the  position  of  r^  is  13.4  mm,  and  the 
position  of  r30  is  14.9  mm.  The  value  of  the  overpressure  at 
14  mm  is  obtained  from  a  plot  of  pressure  versus  position  at 
a  specific  time. 

3.  The  code  proceeds  with  the  calculation  in  the 
following  manner: 

(a)  A  time  step  based  on  the  local  speed  of  sound 
in  the  medium  is  calculated  for  each  cell.  Appendix  C  pro¬ 
vides  a  derivation  of  this  time  step.  The  minimum  of  each 
of  these  time  steps  is  used  as  the  time  for  the  iteration. 

(b)  Acceleration  values  for  each  of  the  cell 
boundaries  are  calculated  and  subsequently  the  velocities  of 
each  cell  boundary  are  calculated  by  multiplying  the 
acceleration  of  the  boundary  by  the  minimum  time  step. 

(c)  New  cell  boundary  position  values  are  calcu¬ 
lated  by  multiplying  the  velocity  values  by  the  minimum  time 
step. 

(d)  New  cell  pressure,  internal  energy,  and  arti¬ 
ficial  viscosity  values  are  calculated  based  on  the  Griineisen 
equation  of  state. 

(e)  The  results  of  that  particular  iteration  are 
printed  out  in  the  form  illustrated  in  Appendix  E. 

4.  The  code  at  present  proceeds  through  2000  iterations. 
Hard  copy  output  is  generated  each  ten  iterations  in  the 
current  computer  system  to  avoid  code  failure  due  to  ex¬ 
cession  of  a  5000  line  output  line  limit.  If  an  alternate 
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system  using  EYESHOC  has  no  such  limit,  hard  copy  output  can 
easily  be  obtained  for  each  iteration. 

5.  From  the  output  generated  by  EYESHOC,  to  obtain  the 
overpressure  at  the  retina,  the  pressure  p  is  plotted  versus 
the  position  r,  and  a  value  of  the  pressure  at  14  mm  is 
obtained.  From  this  value,  the  ambient  pressure  pQ  is  sub¬ 
tracted,  which  leaves  a  value  for  the  overpressure.  This 
value  is  converted  from  ergs  per  cubic  centimeter  (erg/cm3) 
to  atmospheres  (atm)  by  using  the  conversion  factor 

1  atm  *  1.01  x  106  erg/cm3  (22) 

Severu 1  examples  of  pressure  versus  position  plots  at 
various  times  are  illustrated  in  Figure  8.  Plots  A,  B,  and 
C  reflect  the  shock  wave  nature  initially  present  in  the 
pressure  pulse  as  it  propagates.  Plot  B  shows  the  process 
of  "shocking-up"  in  which  the  line  between  the  undisturbed 
gas  and  the  shock  wave  is  essentially  vertical.  As  the  wave 
progresses  to  plot  C,  it  is  beginning  to  relax  from  a  true 
shock  wave  and  degenerate  into  an  acoustic  pulse  propagating 
at  the  local  speed  of  sound.  Plots  D,  E,  and  F  also  serve 
to  illustrate  this  degeneration  into  an  acoustic  pulse.  In 
essence,  EYESHOC  predicts  that  while  a  shock  wave  is  present 
early  in  the  propagation  of  the  wave,  prior  to  reaching  the 
retina  it  has  degenerated  into  an  acoustic  wave  with  the 
pressure  of  the  wave  greatly  reduced  in  magnitude. 

Initial  Condi t ions 

In  order  to  use  EYESHOC  to  predict  overpressure 


rg/cw 


Figure  8.  Fressure  Versus  Distance  for  Various  Times 


levels  at  the  retina,  only  the  specific  internal  energy  and 
pressure  values  for  the  first  cell  need  to  be  calculated  and 
inserted  into  the  code  by  the  user.  Appendix  D  provides  a 
sample  calculation  of  these  values,  and  Table  I  lists  the 
first  cell  specific  internal  energy  and  pressure  values  which 
were  used  in  the  three  test  runs  conducted  with  the  code. 

While  these  values  are  a  function  of  cell  size,  it  was  assumed 
that  the  main  variable  of  interest  was  the  pulse  energy  and 
that  this  energy  was  spread  uniformly  over  the  entire  first 
cell  so  that  in  essence,  the  cell  size  is  not  of  consequence. 
This  dependence  on  pulse  energy  is  also  seen  i,n  the  treatment 


of  Sedov  (Ref  25:238). 


TABLE  I 


EYESHOC  First  Cell  Specific  Internal  Energy 
and  Pressure  Values  for  Test  Runs 


Pulse  Energy  (mJ) 


Internal  Energy  . 
(erg/gm) 


Pressure  (erg/cm3) 


3.75 


7.50 


8.89  x  109 
1 .78  x  10l 0 


1.97  x  10® 
3.93  x  10® 


22.50 


5.33  x  1010 


1.18  x  109 


Code  Val idat ion 

In  an  effort  to  accurately  establish  the  validity  of 
EYESHOC,  the  results  obtained  for  an  ideal  gas  with  a  value 
of  y  equal  to  1.4  were  compared  with  the  results  obtained  by 
Sedov  (Ref  25:247-248)  for  the  identical  case.  The  results 


of  this  investigation  are  discussed  in  detail  and  presented 
in  Appendix  G. 

Figure  9  shows  the  relationship  between  two  dimension¬ 
less  ratios  defined  by  Sedov.  The  ratio  p2/p^  refers  to  the 
ratio  of  the  maximum  pressure  of  the  shock  to  the  ambient 
pressure  while  the  ration  r2/r°  refers  to  the  ratio  of  the 
radial  distance  at  the  point  of  maximum  pressure  to  a  con¬ 
stant  dependent  upon  the  initial  laser  pulse  energy  and 
the  ambient  pressure  of  the  medium.  For  the  values  used  in 
this  thesis,  r°  varied  from  0.33  cm  to  0.60  cm  and  the  ratio 
of  r2/r°  for  the  experimental  distances  to  be  discussed  in 
Chapters  V  and  VI  varied  from  1.1  to  4.3.  For  those  values, 
the  plot  of  versus  T^/r°  predicts  a  radial  dependence 

of  l/r° ' 5  for  y  *  1.022  and  1/r0*06  for  y  »  1.4.  However, 
it  must  be  noted  that  this  is  an  order  of  magnitude  change 
in  the  radial  dependence  and  is  entirely  due  to  the  y  scaling 
present  in  the  theory  of  shock  waves.  If  the  value  of  y  for 
water  calculated  for  this  thesis  is  slightly  wrong,  it  could 
change  the  plot  in  Figure  9  by  a  wide  margin.  The  most 
important  point  to  note  is  that  the  radial  dependence  is 
decreasing  as  the  pressure  drop  becomes  less  and  less 
significant.  This  is  indicative  of  the  great  drop  in 
pressure  as  the  shock  wave  degenerates  into  an  acoustic 
pulse  and  loses  much  of  its  initial  radial  dependence. 

Predicted  Results 

As  shown  above,  three  test  runs  were  made  using  the 


EYESHOC  code,  with  energies  of  3.75  mJ,  7.50  mJ,  and  22.50 
mJ  at  the  focal  point.  The  overpressure  level  at  a  distance 
of  14  mm  from  the  point  of  laser  energy  deposition  was 
determined  in  the  manner  previously  discussed,  and  the  pre¬ 
dicted  values  of  the  overpressures  at  the  retina  are  listed 
in  Table  II. 

TABLE  II 

Predicted  Overpressures  at  the  Retina  for 
EYESHOC  Test  Runs 


Pulse  Energy  (mJ) 

Overpressure  (x  10“ 2  atm) 

3.75 

1.8 

7.50 

3.6 

22.50 

11.0 

We  now  have  overpressures  predicted  based  on  the 
theoretical  modeling  of  a  shock  wave  in  the  eye.  While  the 
actual  measurement  of  overpressures  in  an  eye  is  beyond  the 
scope  of  this  thesis,  we  will  now  discuss  the  test  system 
used  to  model  an  ocular  system,  and  the  overpressures 
measured  by  this  experimental  arrangement. 


V.  Experimentation 


At  this  point,  the  main  thrust  of  this  thesis  moves 
from  theoretical  predictions  to  experimental  investigations. 
This  chapter  will  present  the  materials  employed  and  the 
procedures  followed  in  an  effort  to  obtain  measured  over¬ 
pressure  levels  caused  by  laser  induced  shock  waves.  With 
the  values  of  overpressure  levels  listed  in  Table  II  and 
those  reported  in  Chapter  III,  the  need  was  great  to  attempt 
to  validate  the  predictions  of  EYESHOC  experimentally.  The 
following  describes  how  this  was  accomplished. 

System  Design 

The  discussion  of  the  experimental  system  used  in  these 
experiments  will  be  subdivided  into  discussions  of  the  laser 
system,  the  pressure  transducer  and  associated  electronics, 
and  the  test  cell  and  vitreous  model.  Figure  10  provides 
a  schematic  drawing  of  the  overall  system  while  Figure  11  is 
an  actual  photograph  of  the  system  during  its  use. 

Laser  System 

The  laser  system  used  in  this  study  was  an  entirely  self 
contained  Q-switched  or  mode-locked  Nd: YAG  laser,  assembled 
by  Quantel  International.  The  laser  is  water-cooled,  and 
emits  monochromatic  pulses  of  wavelength  equal  to  1064  nm. 

The  laser  radiation  is  coupled  to  an  ophthalmic  slit  lamp 
via  an  articulated  arm,  which  was  assembled  by  American 


Figure  10.  Test  System  Schematic 
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Figure  11.  Test  System  Photograph 

Medical  Optics,  Incorporated,  so  that  the  output  of  the 
overall  laser  system  is  at  the  rear  of  the  slit  lamp  binocu¬ 
lars,  and  can  be  focused  into  a  patient's  eye.  In  addition 
to  the  Nd:YAG  radiation  exiting  the  system,  the  laser  houses 
a  helium-neon  (He-Ne)  alignment  laser  and  a  visible  light 
beam.  The  internal  optics  are  adjusted  so  that  all  three 
beams  can  be  focused  simultaneously,  i.e.,  if  the  visible 
beam  and  the  He-Ne  alignment  beam  are  in  focus,  the  point  of 
their  focus  is  the  point  of  focus  of  the  Nd:YAG  laser  also. 
The  slit  lamp  binoculars  feature  a  filter  for  both  the  He-Ne 
and  Nd:YAG  beams  to  prevent  injury  to  the  ophthalmologist 
during  the  laser  assisted  psoterior  capsulotomy. 

The  Q-switch  apparatus  of  the  laser  will  be  examined  in 


detail.  However,  the  mode-locking  mechanism  of  the  laser 


system  will  not  be  discussed,  since  no  experiments  were  con¬ 
ducted  for  this  thesis  using  the  laser  system  in  the  mode- 
locked  operating  mode.  Figure  12  provides  a  schematic  of  the 
interior  of  the  laser  cavity  itself,  with  a  list  of  the  optics 
contained  therein. 

In  the  Q-switched  mode  the  laser  emits  15  nsec  pulses 
of  varying  energy  levels  ranging  from  1  mJ  to  35  mJ.  The 
beam  is  near  diffraction  limited  with  a  0.7  mrad  half-angle 
beam  divergence.  The  cavity  is  formed  between  a  flat  100% 
reflecting  mirror  and  a  convex  100X  reflecting  mirror  with 
radius  of  curvature  of  3  m.  These  are  items  1  and  9  respec¬ 
tively  in  Figure  12.  To  insure  Q-switched  operation,  the 
turning  prism,  item  5,  must  be  inserted  into  the  beam  path. 

The  free  running  mode  of  the  laser  in  suppressed  by  the 
intra  cavity  quarter-wave  plate,  WP-1 .  This  element  ensures 
a  90°  rotation  of  the  beam  on  a  double  pass  and  subsequent 
rejection  of  this  mode  by  the  adjacent  Gian -Taylor  polarizer. 
Lasing  in  the  Q-switched  mode  occurs  when  a  quarter-wave  vol¬ 
tage  is  impressed  across  the  Pockels  Cell,  item  2,  which 
further  rotates  the  beam  by  90°.  This  restores  the  beam  to 
the  horizontally  polarized  state,  and  it  propagates  to  the 
laser  rod  for  amplification.  The  beam  undergoes  a  double 
pass  through  the  quarter-wave  plate  WP-2,  and  upon  returning 
to  the  second  Glan-Taylor  polarizer,  the  vertical  component 
of  the  beam  is  diverted  and  used  as  the  laser  output.  The 
horizontal  component  of  the  beam  continues  to  the  Q-switch 
leg  of  the  path  to  sustain  oscillation  in  the  cavity 


(Ref  30:3-5). 

The  laser  was  operated  via  a  footswitch  which  enabled 
the  user  to  fire  the  laser  at  will.  In  addition,  the  laser 
has  the  capability  of  operating  at  a  repetition  rate  of 
10  Hz,  for  te^  purposes.  When  testing  the  laser,  a  turning 
mirror  is  brought  into  the  beam  so  that  the  pulses  can  be 
focused  into  a  Quantel  Energy  Meter  head  and  the  laser  output 
can  be  read  on  an  analog  scale  permanently  attached  to  the 
laser  system  itself. 

This  laser  system  was  the  prototype  for  a  marketable 
system  currently  available  from  American  Medical  Optics,  Inc. 
The  laser  was  made  available  for  use  in  this  study  by 
Dr.  Richard  H.  Keates,  Professor  of  Ophthalmology,  of  the 
Ohio  State  University  Eye  Clinic.  All  tests  were  performed 
at  the  Ohio  State  University. 

Pressure  Transducer  and  Associated  Electronics 

The  pressure  transducer  used  in  these  experiments  was  a 
Kulite  Semiconductor  Pressure  Transducer,  Model  Number  VQH- 
250-5.  A  photograph  of  the  transducer  is  included  in  Figure 
13.  The  transducer  was  a  gage  type,  consisting  of  a  Wheat¬ 
stone  bridge  circuit  with  an  input  impedance  of  2971  ohms  and 
an  output  impedance  of  2505  ohms.  The  transducer  was  rated 
for  a  pressure  of  5  pounds  per  square  inch  (psi)  (14.7  psi  * 

1  atm),  with  a  maximum  pressure  of  20  psi. 

Physically,  the  transducer  has  an  outside  diameter  of 
6  mm  and  is  63  mm  long  overall.  It  consists  of  a  barrel  17  mm 
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Figure  12.  Laser  Oscillator  Schematic  (Ref  30:5) 


long  and  a  nozzle  attached  to  the  rear  of  the  barrel  46  mm 
in  length.  A  diaphragm  is  located  inside  the  barrel  of  the 
transducer.  This  diaphragm  undergoes  a  translation  when  a 
pressure  pulse  is  applied  to  the  transducer.  The  movement 
of  this  diaphragm  then  alters  the  resistance  in  one  leg  of 
the  bridge  circuit,  which  is  eventually  read  out  as  the 
output  of  the  transducer.  The  transducer  was  calibrated  to 
a  sensitivity  of  11  mV/psi,  so  that  an  overpressure  of 
magnitude  equal  to  1  psi  reaching  the  diaphragm  would  cause 
an  11  mV  signal  to  be  generated  by  the  bridge  circuit.  Thii 
signal  could  then  be  detected  by  an  oscilloscope,  or  as  was 
used  in  this  study,  by  a  digital  or  analog  voltmeter. 


Figure  13.  Pressure  Transducer 


The  power  supply  for  the  transducer  consisted  of  an 
Eveready  BA-53  dry  cell  in  series  with  a  500  ohm  variable 


resistor.  This  type  of  voltage  source  was  used  in  order  that 
the  20  V  DC  excitation  voltage  which  was  used  as  the  calibra¬ 
tion  parameter  for  the  transducer  could  be  maintained  at  a 
stable  level.  In  addition ,  the  use  of  a  dry  cell  as  a  voltage 
source  decreases  the  overall  noise  level  of  the  system,  which 
would  have  been  greater  had  the  voltage  source  been  a  conven¬ 
tional  line  DC  power  supply.  The  dry  cell  and  variable 
resistor  were  housed  in  an  aluminum  shield  to  further  lower 
the  noise  level,  and  the  voltage  from  the  dry  cell  to  the 
transducer  was  delivered  via  shielded  cable  using  BNC  type 
connections. 

The  output  of  the  transducer  was  measured  using  a  Data- 
Precision  3500  Digital  Voltmeter.  The  voltmeter  was  used  in 
a  mode  such  that  the  minimum  detectable  voltage  was  0.1  mV. 

The  transducer  was  also  housed  in  an  aluminum  shield,  which 
was  affixed  to  the  rear  of  the  test  cell,  with  the  barrel  of 
the  transducer  entering  the  test  cell  through  a  hole  bored  in 
the  back  wall  of  the  cell.  The  output  of  the  transducer  was 
connected  to  the  voltmeter  again  via  shielded  cable  using  BMC 
type  connections. 

Test  Cel 1  and  Vitreous  Model 

The  experiments  conducted  in  this  study  employed  a  glass 
and  plexiglas  test  cell  developed  by  J.  Riggins  (Ref  5:57). 

The  cell  is  illustrated  in  Figure  14.  The  cell  was  constructed 
of  plexiglas  on  the  sides,  top,  and  bottom,  while  tne  two  faces 


—  front  und  back  --  in  the  path  of  the  beam  were  constructed 
of  glass.  The  top  of  the  cell  was  removable  to  allow  the  cell 
to  be  filled  with  various  substances  for  testing  purposes. 

The  lid  was  also  slotted  to  allow  for  the  translation  of  a 
small  plexiglas  cylinder,  which  was  used  to  aid  in  the  assign¬ 
ment  of  distances  from  the  transducer.  The  slide  which  moved 
the  cylinder  was  kept  under  tension  by  two  springs  while 
movement  was  obtained  by  a  fine  threaded  screw  which  provided 
movement  of  0.625  nun  per  rotation. 


Figure  14.  Test  Cell 

For  this  study,  the  test  cell  was  filled  with  balanced  salt 
solution  (BSS).  This  solution  had  previously  been  shown 

55 


nXMr 


.. nA  v.vt V-. %  N.;. _v_ *.  . . v 


to  have  approximately  the  same  breakdown  characteristics  of 
the  vitreous  itself  (Ref  5:85).  Actual  vitreous  was  originally 
to  have  been  employed  for  testing  purposes,  however  human 
vitreous  in  the  quantities  needed  —  50  ml  per  test  —  was 
virtually  unobtainable.  The  use  of  hog  vitreous  was  consid¬ 
ered,  however  the  material  obtained  was  not  of  sufficient 
purity  to  lend  credence  to  any  measurements  which  could  have 
been  generated. 

The  overall  dimensions  of  the  test  cell  were  40  mm  high 
by  40  mm  wide  by  32  mm  deep.  These  dimensions  yield  a  volume 
of  51.2  ml  for  the  test  cell.  As  was  mentioned,  the  front 
and  back  of  the  cell  were  constructed  of  glass.  The  glass 
front  of  the  cell  was  removed  in  order  to  determine  the 
absorption  of  the  laser  radiation  due  to  the  glass  itself. 

The  glass  was  placed  in  a  Perk in-Elmer  Model  350  Spectro¬ 
photometer,  and  a  plot  of  transmittance  versus  wavelength 
from  0.6  to  2.7  microns  was  generated.  This  plot  is  shown 
in  Figure  15.  To  identify  the  NdrYAG  laser  wavelength 

absolutely,  a  Laser  Optics  narrow-bandpass  filter  selectively 

o 

coated  with  a  100  angstrom  (A)  coating  for  1.06  microns  was 
used.  The  results  of  this  investigation  revealed  that  75X  of 
the  incident  radiation  at  1.06  microns  was  being  transmitted 
to  the  interior  of  the  cell.  Using  this  data,  all  laser 
pulse  energies  were  corrected  for  this  factor  in  subsequent 
uses. 


The  laser  system,  the  pressure  transducer,  and  the  test 
cell  filled  with  balanced  salt  solution  were  integrated . into 
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Figure  15.  Cell  Front  Transmission  Data 


a  total  test  system  for  the  determination  of  the  laser 
induced  shock  wave  overpressures  at  specific  distances  from 
the  diaphragm  of  the  transducer.  Figure  16  is  a  close-up  of 
the  test  cell  with  the  transducer  attached  to  the  rear  of  the 
cell.  The  entire  cell-transducer  assembly  is  attached  to  the 
output  end  of  the  laser  system  for  actual  experimentation. 

Ex per imental  Procedures 

Once  the  equipment  was  assembled,  the  following  steps 
were  performed  in  order  to  determine  the  overpressures  caused 
by  several  laser  pulse  energies  at  varying  distances  from  the 
pressure  transducer: 


Figure  16.  Test  Cell  with  Attached  Pressure  Transducer 


1.  The  output  of  the  laser  was  diverted  into  the 
energy  meter  head,  and  the  pulse  energy  was  set  to  a  value 

of  5  mJ.  Subsequent  tests  were  conducted  with  pulse  energies 
of  10  mJ,  and  30  mJ,  as  recorded  using  the  readout  of  the 
energy  meter.  Making  use  of  the  75%  transmission  factor, 
these  values  correspond  to  pulse  energies  at  the  focal  point 
of  3.75  mJ,  7.50  mJ,  and  22.50  mJ,  respectively. 

2.  The  output  of  the  laser  was  again  allowed  to  enter 
the  test  cell.  The  translating  cylinder  was  initially  set 
at  a  distance  of  14  mm  from  the  transducer,  but  was  moved 

to  distances  of  10  mm,  8  mm,  and  6  mm  from  the  transducer  in 
subsequent  experiments. 

3.  Ten  laser  pulses  were  focused  at  the  distances 
mentioned  above,  and  the  overpressure  induced  by  each  was 
measured  using  the  digital  voltmeter.  The  average  over¬ 
pressure  for  each  energy  and  distauce  combination  was 
determined  by  calculating  the  mean  of  the  ten  values,  and 
the  error  in  the  measurements  was  calculated  by  determining 
the  standard  deviation  of  each  of  the  mean  values.  The  dis¬ 
tances  were  determined  to  within  +  1  mm  using  the  translating 
cylinder,  the  slit  lamp  binoculars,  and  the  spark  which  is 
visible  upon  breakdown  of  the  BSS.  It  should  be  noted  that 
each  of  the  ten  pulses  focused  into  the  test  cell  initiated 
breakdown  of  the  BSS. 

4.  After  each  energy  and  distance  combination  had  been 
tested,  the  test  cell  was  re-filled  with  fresh  BSS  to  prevent 
the  accumulation  of  any  possible  artifacts  created  in  the 
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breakdown  process.  This  was  done  to  maintain  as  much 
continuity  in  the  measurements  as  possible. 

This  concludes  the  discussion  of  the  test  system,  and 
the  manner  in  which  the  tests  were  conducted.  The  next 
chapter  provides  the  results  obtained  in  these  experiments 
and  compares  these  results  with  those  predicted  by  theory. 


VI.  Results  and  Discussion 


This  chapter  will  present  the  overpressure  measured  in 
the  experiments  described  in  Chapter  V  and  compare  these 
values  to  the  predictions  given  in  Chapter  IV.  In  addition, 
the  results  will  be  examined  in  order  to  determine  whether 
or  not  they  lit  the  model  of  a  spherical  shock  wave  with 
reference  to  pressure  versus  energy  and  pressure  versus 
distance. 

Overpressure  versus  Energy 

The  overpressures  measured  at  the  transducer  for  varying 
distances  from  the  transducer  as  a  function  of  input  pulse 
energy  are  listed  in  Table  III.  It  can  be  seen  that  over¬ 
pressures  ranged  from  0.018  atm  to  0.432  atm,  or  changes  in 
the  overall  pressure  at  the  transducer  of  1.7%  and  42. 1%  re¬ 
spectively.  If  the  values  obtained  for  a  distance  of  14  mm 
are  examined,  it  can  be  seen  that  the  maximum  overpressure 
was  0.077  atm,  or  an  overall  change  of  7.5%  from  ambient 
pressure.  If  the  BSS  test  cell  environment  is  a  valid  sub¬ 
stitute  for  the  eye,  which  we  believe  it  to  be,  then  these 
would  be  the  overpressures  expected  during  the  posterior 
capsuLotomy  procedure.  It  is  very  important  to  note  that 
the  largest  input  pulse  energy  did  not  alter  the  pressure 
by  more  than  10%  at  this  distance. 

Based  upon  theoretical  treatments  such  as  by  Sedov 


the  overpressure  should  be  linear  when  plotted  versus  the 


input  pulse  energy.  Figure  17  is  a  plot  of  the  values 
listed  in  Table  III,  for  each  distance  from  the  transducer. 

The  correlation  coefficient  was  calculated  for  each  of  the 
lines,  with  the  result  being  that  each  of  the  four  lines 
plotted  in  Figure  17  had  a  correlation  of  0.999  or  greater. 

A  perfect  correlation  between  data  points  would  yield  a 
correlation  coefficient  equal  to  one,  hence  it  is  obvious 
that  the  overpressures  measured  in  this  study  are  indeed 
linear  with  the  pulse  energy,  as  expected. 

Overpressure  versus  Distance 

Based  on  the  degeneration  of  the  shock  wave  into  an 
acoustic  pulse  the  expected  variation  of  the  overpressure 
with  distance  from  the  laser  induced  breakdown  site  to  the 
transducer  is  for  the  overpressures  to  show  a  1/r2  fall-off 
as  the  distance  increases.  In  view  of  this,  if  the  over¬ 
pressures  are  plotted  versus  the  distance,  to  accurately  show 
this  type  of  trend,  we  should  multiply  the  overpressures  by 
the  square  of  the  distances. 

The  values  of  the  overpressures  obtained  versus  the  dis¬ 
tance  from  the  transducer  are  listed  in  Table  IV.  In  addition 
these  values  are  plotted  in  Figure  18,  which  gives  three  plots 
of  overpressure  times  distance  squared  versus  distance,  one 
plot  for  each  of  the  three  energies  tested.  As  is  noted  on 
the  figure,  the  broken  line  indicates  the  measured  values, 
while  the  solid  line  indicates  the  results  of  a  least-squares 
fit  analysis  of  the  data.  As  can  be  seen,  the  values  obtained 


Figure  17.  Overpressure  versus  Energy 


in  the  least-squares  analysis  are  virtually  linear,  as  is 
expected  based  on  the  increase  in  energy.  That  is,  the  value 
for  E  *  22.75  mJ  is  approximately  six  times  the  value  for 
E  =  3.75  mJ.  From  the  graphical  plots  in  Figure  18,  it 
seems  indeed  safe  to  conclude  that  the  pressure  does  show  a 
1/r2  fall-off  with  distance. 


Comparison  of  Predicted  versus  Measured  Overpressures 

The  magnitudes  of  the  measured  and  predicted  over¬ 
pressures,  and  the  percent  change  from  ambient  pressure  which 
they  represent  are  listed  in  Table  V.  The  differences  in 
measured  and  predicted  variations  from  ambient  pressures  can 
be  seen  to  be  less  than  5%  in  each  case  —  assuming  the 
worst  case  analysis  in  which  the  lower  range  of  the  standard 
deviation  is  used  as  the  value  for  the  measured  overpressures 
These  values  represent  the  maximum  pressure  levels  present 
in  an  acoustic  pulse  propagating  at  the  local  speed  of  sound. 
It  should  be  noted  that  these  values  are  for  the  breakdown- 
to-transducer  distance  of  14  mm  only.  The  values  listed  in 
Table  V  are  plotted  in  Figure  19. 

If  the  trends  exhibited  in  Figure  19  are  analyzed 
closely,  it  can  be  seen  that  the  difference  between  the  two 
lines  increases  as  the  energy  increases.  This  is  an  indica¬ 
tion  of  a  discrepancy  occuring  in  either  the  code,  the 
measurements,  or  possibly  in  both.  If  the  plot  is  examined 
closely,  it  can  be  seen  that  the  line  indicating  the  pre¬ 
dicted  overpressures  is  completely  linear,  while  the  line 


indicating  the  measured  overpressures  begins  to  show  a  fall- 
off  from  its  original  slope,  therby  deviating  from  total 
linearity  as  the  energy  of  the  pulse  increases.  However,  if 
the  deviation  of  the  measured  values  is  considered,  it  can  be 
seen  that  for  E  =  7.5  mJ,  the  predicted  value  lies  in  the 
range  of  the  measured  v?.lue  and  for  E  *  22.50  mJ,  the  differ¬ 
ence  between  the  largest  measurement  and  the  predicted  value 
is  only  0.011  atm,  which  is  a  miniscule  difference.  At  the 
levels  of  overpressures  being  discussed  in  this  thesis,  such 
differences  cannot  be  considered  significant.  The  measured 
values  for  this  case  are  deemed  quite  good,  in  view  of  the 
fact  that  at  the  greatest  pulse  energy,  the  greatest  deviation 
between  the  values  yields  a  difference  of  4.4%  in  the  percent 
change  from  ambient. 

This  concludes  the  presentation  and  discussion  of  the 
results  obtained  in  the  experimental  portion  of  this  investi¬ 
gation.  At  this  point,  hwe  have  examined  the  theory  associated 
with  both  the  ocular  system  in  humans  and  the  physics  of  shock 
waves;  we  have  derived  a  theoretical  hydrodynamic  code  for 
predicting  laser  induced  shock  wave  overpressures  at  the 
retina;  and  we  have  developed  and  implemented  a  test  system 
designed  to  measure  such  overpressures.  All  that  remains  is 
to  assemble  the  data  generated  from  both  the  analytical  and 
experimental  portions  of  this  study  and  draw  conclusions  from 
it.  This  is  the  topic  of  Chapter  VII. 


VII.  Recommendations  and  Conclusions 

The  purpose  of  this  chapter  is  to  analyze  the  findings 
in  this  study  from  an  overall  point  of  view.  We  wish  to 
examine  what  the  study  has  accomplished,  and  recommend 
further  investigations  in  the  area  of  Nd : YAG  laser-assisted 
ocular  surgery.  In  addition,  clinical  conclusions  will  be 
presented  in  this  chapter  which  can  be  used  as  possible 
guidelines  for  the  ophthalmologist  when  performing  posterior 
capsulotomies  with  Nd:YAG  lasers. 

This  study  examined  the  shock  wave  generated  in  the 
posterior  chamber  of  the  human  eye  during  posterior  capsulo¬ 
tomies  performed  with  Nd:YAG  lasers.  These  shock  waves  were 
modeled  using  a  one-dimensional  hydrodynamic  code  in  spherical 
coordinates  utilizing  a  Lagrangian  mesh.  The  code  predicted 
that  the  shock  waves  initially  present  degenerated  into  acous¬ 
tic  pulses  prior  to  reaching  the  retina  and  that  the  over¬ 
pressures  at  the  retina  are  less  than  0.15  atm.  In  addition, 
overpressures  caused  by  Nd:YAG  laser  systems  currently  used 
by  several  ophthalmologists  to  treat  patients  were  measured 
in  a  test  chamber  containing  balanced  salt  solution.  These 
measurements  were  made  for  varying  distances  from  the  point 
of  energy  deposition  to  the  pressure  transducers,  as  well  as 
for  varying  input  laser  pulse  energies.  The  overpressures 
measured  in  this  experiment  varied  from  0.018  atm  to  0.432 
atm.  A  comparison  of  predicted  overpressure  values  obtained 
from  the  code  and  measured  overpressure  values  obtained  from 


the  experiments  indicated  a  high  degree  of  agreement  between 
the  two  methods  of  generating  the  overpressure  at  the  retina 
during  this  type  of  ocular  surgery. 

Further  studies  and  experiments  are  necessary  for  a 
complete  understanding  of  the  processes  involved  when  the 
energy  from  a  Nd:YAG  laser  is  focused  into  the  posterior 
chamber  of  the  eye.  Such  studies  should  include  both  addi¬ 
tional  experimental  and  theoretical  investigations.  To  this 
end,  changes  are  proposed  for  the  shock  code  and  future 
experiments  are  suggested. 

Shock  Code  Limi tat  ions  and  Proposed  Changes 

1.  This  thesis  was  concerned  exclusively  with  the  retina 
in  terms  of  the  effects  of  the  shock  wave.  Analytical  pre¬ 
dictions  such  as  those  obtained  from  EYESHOC  should  be  extended 
to  other  possibly  sensitive  portions  of  the  eye  to  determine 
the  overpressures  at  these  tissues.  Examples  include  the  iris, 
cornea  —  assuming  no  lens  present  —  or  the  lens  itself. 

This  would  be  a  simple  calculation  based  on  the  results  presented 
in  this  study  if  the  distances  from  the  point  of  laser  energy 
deposition  to  the  tissues  of  interest  were  known.  These  cal¬ 
culations  should  be  made,  and  subsequent  codes  devised  to 
determine  their  accuracy. 

2.  Unfortunately,  EYESHOC  suffered  from  some  inherent 
numerical  instabilities  which  survived  all  efforts  to  remove 
them.  Such  instabilities  are  not  totally  uncommon  in  codes 


of  this  type  (Ref  29:34-66),  however  with  extensive  work 
they  can  be  reduced  and/or  eliminated  altogether.  The 
primary  aid  in  reducing  instabilities  is  through  the  use  of 
artificial  viscosity,  q,  which  is  designed  to  smooth  out  the 
shock  front  and  prevent  such  instabilities.  More  extensive 
analysis  of  the  artificial  viscosity  present  in  the  code 
should  be  performed. 

3.  Finally,  an  effort  should  be  made  to  incorporate 
the  actual  viscosity  of  the  vitreous  into  the  analytical 
model  designed  for  EYESHOC.  The  vitreous  is  four  times  as 
viscous  as  water  (Ref  9:252),  but  this  does  not  appear  in 
the  model,  which  is  based  almost  exclusively  on  water.  This 
could  possibly  appear  in  relation  to  the  density,  however 
this  author  offers  no  suggestion  as  to  the  manner  of  incor- 
portaiton . 

The  three  points  listed  herein  should  be  the  topics  of 
future  theoretical  analyses.  However,  it  cannot  be  over¬ 
looked  that  the  results  of  EYESHOC  are  in  very  close  agree¬ 
ment  with  the  overpressures  obtained  in  experimental  studies, 
booh  by  other  experimenters  and  those  accompanying  this  thesis 
effort.  In  view  of  this,  it  must  be  concluded  that  EYESHOC 
is  a  useful  tool  for  predicting  the  overpressures  associated 
with  the  Nd:YAG  laser  procedure,  and  in  this  respect,  EYESHOC 
accomplishes  its  intended  purpose. 

Future  Experimentation 


The  experiments  conducted  in  this  thesis,  as  well  as 


those  discussed  in  Chapter  II  lay  a  firm  experimental  found¬ 
ation  needed  to  characterize  NdrYAG  laser  assisted  posterior 
capsulotomies.  However,  future  experiments  are  necessary 
and  should  be  of  two  main  types,  in  vivo  and  in  vitro  studies. 

Initially,  the  next  group  of  experiments  should  be  per- 
formed  in  an  ijn  vitro  manner,  i.e.  with  a  test  cell  and 
vitreous  model  as  was  used  in  this  thesis.  However,  two  main 
changes  should  be  instituted  in  these  studies.  First,  the 
test  cell  should  be  altered  in  geometry  so  that  it  is  the 
same  size  and  shape  as  an  actual  eye.  This  cell  should  be 
fitted  with  multiple  pressure  transducers  rather  than  the 
single  transducer  used  in  this  study.  This  would  enable 
overpressure  measurements  to  be  obtained  for  areas  of  the  eye 
in  addition  to  the  retina.  Second,  a  substance  other  than 
BSS  should  be  used  for  the  vitreous  model.  This  will  involve 
preliminary  experimentation  to  identify  a  suitable  model  for 
the  vitreous  which  not  only  exhibits  the  same  breakdown 
characteristics,  but  is  approximately  of  the  same  viscosity 
as  the  vitreous.  These  two  changes  would  enable  the  closest 
correlation  between  a  test  system  and  an  actual  patient  to  be 
achieved . 

Subsequently,  experiments  should  be  undertaken  iii  vivo 
with  both  rabbits  and  monkeys,  to  measure  the  overpressures 
generated  in  living  systems.  As  mentioned  above  for  test 
cell  studies,  multiple  pressure  trnasducers  should  be  im¬ 
planted  in  these  animals  ,  and  values  of  the  overpressures 


should  be  determined  at  all  possible  points  in  the  eye.  To 
totally  investigate  the  effects  of  such  shock  waves,  electron 
microscopy  should  be  performed  on  all  ocular  tissues  from 
the  animals  irradiated  in  the  manner  discussed  in  this 
thesis. 

Finally,  these  same  studies  should  be  conducted  on 
human  eyes  which  have  been  obtained  from  an  eye  bank.  This 
would  of  course  simulate  the  actual  conditions  present  ' 
during  a  posterior  capsulotomy  as  closely  as  possible. 


Clinical  Conclusions 


The  most  important  goal  of  this  thesis  was  to  develop 
guidelines  for  the  ophthalmologist  to  use  when  performing  a 
posterior  capsulotomy  on  a  patient  with  regard  to  the  over¬ 
pressure  at  the  retina  caused  by  a  laser  induced  shock  wave. 
To  this  end,  we  have  determined  the  values  of  such  over¬ 
pressures,  and  now  define  a  maximum  overpressure  which 
should  not  be  exceeded. 

Based  on  the  results  obtained  in  this  thesis,  a  maximum 
level  of  acceptable  overpressure  of  0.10  atm  is  proposed  as 
the  treatment  threshold  when  performing  posterior  capsuloto- 
mies  with  Nd:YAG  lasers.  This  value  would  allow  laser 
pulses  of  up  to  10  mJ  to  be  used  with  no  ill  effects.  These 
pulses  should  be  of  nanosecond  pulsewidth,  as  were  those 
employed  in  this  study.  In  addition,  this  overpressure 
level  will  allow  for  pulse  energy  deposition  as  far  as  6  mm 


behind  the  posterior  capsule  without  damage  to  the  retina. 
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This  would  enable  the  ophthalmologist  to  utilize  the  Nd:YAG 
laser  in  the  removal  of  vitreous  bands  as  well,  without 
retinal  damage  concerns. 

One  of  the  primary  reasons  for  assigning  the  0.10  atm 
level  as  the  maximum  acceptable  overpressure  deals  with  the 
current  use  of  Nd:YAG  lasers  by  practicing  ophthalmologists. 
Several  thousand  patients  have  undergone  posterior  capsulo- 
tomies  performed  by  doctors  Keates,  Fankhauser,  and  Aron-Rosa, 
and  not  one  patient  has  suffered  retinal  damage  attributable 
to  a  laser  induced  shock  wave  or  acoustic  pulse.  The  lasers 
used  by  these  doctors  all  meet  the  criteria  set  forth  above, 
i.e.,  pulses  of  less  than  10  mJ,  nanosecond  pulsewidths,  and 
distances  to  the  retina  normally  greater  than  10  mm  during 
a  posterior  capsuiotomy  procedure.  With  these  limitations 
met,  an  ophthalmologist  could  perform  this  type  of  procedure 
without  fear  of  patient  injury. 
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APPENDIX  A. 


Calculation  of  ^  for  Water 


One  of  the  most  important  requirements  in  the  initial 
stages  of  the  development  of  EYESHOC  was  to  obtain  a  value 
of  y  —  the  ratio  of  the  specific  heat  at  constant  pressure 
to  the  specific  heat  at  constant  volume  —  for  water,  as  this 
was  the  basis  of  the  model  for  the  vitreous.  This  y  factor 
was  required  in  the  Griineisen  equation  of  state.  The  calcu¬ 
lation  is  as  follows,  beginning  with  the  relationship  between 

C  and  C  (Ref  31:168), 
p  v  * 


Cp  -  Cy  =  VTa 2 / k 


(A-l  ) 


where 


V  is  the  volume  (cm3) 

T  is  the  absolute  temperature  (°K) 
a  is  the  volume  coefficient  of  expansion  (°K-1) 
k  is  the  isothermal  compressibility  (cm3/erg) 


For  this  calculation,  these  variables  assumed  the  values 


(Ref  32:F-5)  , 


=  1  cm3 


=  310.15  K 
=  3.62  x  L0~V°K 


=  4.45  x  I0~ixcm3/erg 


Rearranging  Equation  (A-l),  we  obtain, 

Cy  =  Cp  -  VTa2 /<  ( A-2 ) 

and  for  the  absolute  temperature  given, 

Cp  =  4.1755  x  107  erg/gm-°C 
so  that  upon  solution, 

Cy  =  4.0842  x  107  erg/gm-°C 

Finally,  substituting  these  values  of  Cp  and  Cy  into  Equation 
( 13) ,  we  obtain  , 

Y  =  1.022 

At  this  point,  we  should  note  that  all  of  the  calcula¬ 
tions  performed  in  this  thesis  are  highly  dependent  upon  the 
value  of  y  derived  for  water.  If  the  equations  are  examined, 
it  can  be  seen  that  the  dependence  of  the  Griineisen  equation 
of  state  appears  as  (y-1).  For  the  value  of  y  obtained  above, 
this  is  equal  to  2.2  x  10“2 .  If  the  value  o 2  y  is  incorrect 
by  as  much  as  0.01,  the  value  of  (y-1)  is  changed  to  3.2  x  10~2 
which  is  an  increase  of  45X,  while  the  change  in  the  value  of 

Y  is  a  mere  1%.  Extreme  care  was  taken  in  the  calculation  of 

Y  as  presented  herein,  however,  it  must  be  realized  that  even 
if  a  minor  error  is  present  in  this  calculation,  it  will  be 
magnified  by  the  equations  in  use. 
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APPENDIX  B 


Elimination  of  Specif ic  Internal  Energy  in  Time 
Derivative  of  Pressure 

lly,  the  Griineisen  equation  of  state  is  expressed 
density  rather  than  the  specific  volume  V.  Thus 

p  =ph+(y-1)p(E-Eh)  (B-l 

ph  =  pqC2p(p  -  pQ)/[p  -  i  y(p  -  pq>]2  (B-2 

eh  =  i[C(p  -  pQ)/[p  -  £y(p  -  P0)l  ]2  (B-3 


Solving  for  E  in  the  above  equations, 


numerator  =  [ ( y-1 ) ( p-pQ)p  -  \ y( Y-1 )(P-PQ>2  -  PPD 

+  (y-1 )pQ( P-Pq  > ]  +  [-(y-1  )<P-P0)2  +  2(p-pq)po 

+  4y(Y-1 )(P-Pq)2  -  2(y-1)(P-P0)P0] 

numerator  =  pQl (y-4y )(P-P0>  -PQ] 

numerator  =  Poliy(.p-PQ)  -pQ] 

Therefore,  the  time  derivative  of  the  specific  internal 
is  given  by , 


dE  1  dp  , 

dt  (y-1)  dt 


PoC‘ 


(Y-U(P-P0)P0 


(Y-1)PZ  +  (Y-1)  P - ( Y ~ 1  )  ( P- P _ ) I 


dp 

Mt 


Solving  for  dp/dt, 


/  ,  x  dE 

(Y-DP^ 


■  n  ) 

-  p  ♦  fOo0’ 

§Y(P-P0)P0 

[  p-4y(p-p0)P" 

k 

- 

dp 

dt 


or, 


dt 


p  +  (y-1  )(P+q)-P2P„C2 


4y(p-p0)p0 


Ip-  4y(p-p0)] s 


v-u 


(B-9) 

(B-10) 

(B-ll) 

energy 

(B-12) 

(B-13) 

(B-14 ) 


if  the  artificial  viscosity  and  the  value  of  dp/dt  are  included 


AD-A124  731 
UNCLASSIFIED 


NEODVNIUH:  VAG  LASER  INDUCED  SHOCK  NAVE  OVERPRESSURES  “ 
AT  THE  RETINACUT  AIR  FORCE  INST  OF  TECH 
HRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGINEERING 
J  F  TURNER  DEC  82  AFIT7GEP7PH782D-22  F/G  6718  NL 
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APPENDIX  C 


Calculation  of  Time  Step  Based  on  Local 
Adiabatic  Speed  of  Sound 


The  adiabatic  speed  of  sound  of  a  cell  J  at  a  time  t  can 


be  written 


C  =  (dp/dp )‘ 


(C-l) 


which  can  be  expressed  as 


C  -  [  (dp/dt ) / (dp/dt  l]1 


(C-2) 


Substituting  Equations  (1)  and  (B-14)  into  the  above  express¬ 
ion,  we  obtain. 


p  +  (y-1)(P+Q)  - 


P2P0C  [iY(p-P0)P0l 
[p  ■-  iY<P-P0)]  9 


(C-3) 


If  we  assume  that  C  is  constant  over  the  time  interval  of 
interest,  and  if  the  shock  wave  is  only  allowed  to  cross 
one-fifth  of  a  cell  in  a  time  unit  At,  then  this  At  is  given 


PiP0C2[lyip-Q0)p0T1  )  4 


At  -  0.2(Ax)  p  P  +  <Y-D(P+q>  -  [p_  4Y(p_p  )l3  f  (C-4) 


where , 


Ax  is  the  distance  travelled  by  the  boundary  of  the  cell 


It  should  be  noted  that  If  p  *  pQ,  the  value  of  the 
speed  of  sound  reduces  to 

C  -  (YP/P)*  <C-5) 

which  is  the  ideal  gas  value  of  the  speed  of  sound,  without 
the  artificial  viscosity  included,  and  the  time  At  becomes 

At  =  0.2  (Ax)/(yp/p)^ 
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APPENDIX  D 


Calculation  of  First  Cell  Specific  Internal 


Energy  and  Pressure 


For  this  example,  consider  an  energy  present  in  the 
laser  pulse  of  5  mJ: 


E  =  5  mJ  =  5  x  10“  erg 
o 


(D-l) 


If  we  consider  a  sphere  of  radius  100  microns,  then  the 
volume  of  this  sphere  is  given  by 


V  =  |  trr3 


V  =  4.189  x  10“6  cm3 


(D-2) 


To  determine  the  mass  contained  in  this  sphere,  the  volume  is 


multiplied  by  the  density  pQ,  according  to 


m  =  pQV 


CD-3) 


m  =  4.218  x  10"6  gm 


To  calculate  the  specific  internal  energy,  we  divide  the 
pulse  energy  by  the  mass,  so  that 


E  ■  E  /m 
o 


CD-4) 


■  ■r--«  .<*—• Itr-  •* — •L  *.  •  ••  •.  •  •  •. 


E  »  1.185  x  10l°  erg/g 


If  we  examine  the  Gruneisen  equation  of  state,  we  see  that 

at  t  *  t  ,  p  *  p  .  and  therefore. 
o*o 

PH  «  0  (D-5a ) 

Ejj  *  O  (D-5b  ) 


and 


P  *  P0<Y  “  1)E 


CD-6) 


or, 


p  -  2.625  x  10*  erg/cm3 


APPENDIX  E 

EYESHOC  Code  Listing  and  Sample  Output 

This  appendix  provides  the  complete  code  listing  for 
EYESHOC,  along  with  a  complete  list  of  the  variables  used  in 
the  code.  In  addition,  the  output  generated  by  the  program 
is  illustrated  in  Figure  E-l. 


;his  code  models  the  one-oi mens ional  prop AGXTF9N  OF  rmerwm 

IH  THE  HUNAN  EVE*  THE  SHOCK  IS  CAUSED  BV  FOCUSING  A  NOSYAG  LASER 

in  Of  near  the  posterior  capsule  of  the  eve  to  remove  opacities  in 

(HE  CAPSULE*  THIS  COOE  USES  A  LAGRANGI  AN  MESH  TO  DESCRIBE  THE 
PROPAGATION.  THE  SHOCK  IS  ALL  3WEO  TO  PROPAGATE  TO  THE  RETINA* 

WITH  EMPHASIS  BEING  PLACED  ON  THE  AMPLITUDE  OF  THE  PRESSURE  FRONT 
AT  THE  RETINA  IN  AN  EFFORT  TO  QUANTITATIVELY  ESTIMATE  THE  HAXtNUM 
LASER  FLUX  WHICH  CAN  SAFELT  BE  DEPOSITED  IN  THE  POSTERIOR  CAPSULE 
AND  *.OT  CAUSE  RETINAL  DETACHMENT. 

PROGRAM  EvESHOC  (INPUT*  OUTPUT*  TAPE6S0UTPUT) 

DEFINITION  OF  VARIABLES 

A  -  FACTOR  EQUAL  TO  APPROXIMATELY  ONE-HALF  THE  NUMBER  OF  CELLS 
OVER  WHICH  THE  SHOCK  IS  TO  BE  SPREAD 
ALPHA  -  FACTOR  WHICH  DETERMINES  THE  COOE  GEOMETRY*  WHERE 
ALPHA  s  1  FOR  CARTESIAN  COORDINATES 
ALPHA  s  2  FOR  CYLINDRICAL  COORDINATES 
ALPHA  =  i  FOR  SPHERICAL  COORDINATES 
OELT  -  TIME  STEP  FOR  EACH  CYCLE 

OT  -  INDIVIDUAL  TIME  STEP  VALUES*  THE  MINIMUM  OF  WHICH  IS  OELT 
OX  -  WIDTH  OF  EACH  CELL  (CM)  v  * 

t  -  INTERNAL  ENERGY  ( ERG/6M) 

EH  -  INTERNAL  ENERGY  CORRECTION  FACTOR  DEPENOCNT  ON  CHANGE  IN 
OENSITY  DUE  TO  SHOCK  HAVE 

F  -  COMPRESSION  FACTOR  INCLUDED  TO  CONSIDER  BOUNDARY  SHIFTS 
FPRM  -  STORAGE  FILE  FOR  F  VALUE  OF  EACH  PASS 

FRACT  -  EXPERIMENTALLY  DETERMINED  CONSTANT  RANGING  FROM  0.1  TO  O.M 
GAMMA  -  hATIO  OF  SPECIFIC  HEATS  CCP/CV) 

GAMP- M  -  GAMMA  MINUS  ONE 
J  -  CELL  NUMBEP 

OLOV  -  VELOCITY  VALUES  DETERMINED  ON  PREVIOUS  PASS  (CM/SEC) 
uLUX  -  POSITION  VALUES  DETERMINED  ON  PREVIOUS  PASS  (CN> 

P  -  PRESSURE  IN  CELL  J  (E RG/CN*»3) 

t'H  -  PRESSURE  CORRECTION  FACTOR  DEPENDENT  ON  CHANGE  IN  DENSITY 
DUE  TO  SHOCK  NAVE 

•J  -  ARTIFICIAL  VISCOSITY  FACTOR  USED  TO  DEFINE  A  CLEAN  SHOCK  FRONT 
**0  -  DENSITY  OF  CELL  J  (GM/CM**)) 

T  -  TEMPERATURE  OF  CELL  J  (DEG  K> 

TIME  -  TOTAL  ELAPSED  TIME  SINCE  LASER  ENERGY  DEPOSITION  (SEC) 

UOOT  -  ACCELERATION  OF  CELL  J  (CM/SEC**2) 

V  -  VELOCITY  OF  BOUNDARY  OF  CE.L  J  (CM/SEC) 

K  -  BOUNDARY  GF  CELL  J  (CM)  t 

DIMENSION  X(S0)«  P(S0)«  Q(53>*  tf(50)*  R0(S0)«  DT(S0I*  T(S0»*  E(SO) 
DIMENSION  UOOT (SO ) *  OLOX(SO)*  3L0V(S0)*  EH(SO)*  PH(SO)*  OX(SO) 
DIMENSION  Af.UM(50 ) *  uPE£D(50) 

INTEGER  CYCLE 

initialization  of  variables 

A  s  t* 


*  v* 


ALPHA  si. 

FA AC '  s  0.2 
GAMP  M  =  2.2E-2 
*  t ME  s  0. 

OXCl)  S  l.OF-2 
03  5  J  :  2i  50 
5  OXCJ)  s  1.10*OXCJ-1) 

XCl)  S  o. 

DO  10  J  :  2t  45 
DM  J—  1 »  s  0. 

ECJ-1)  S  4.64C7 
PCJ-1)  s  1.0AE6 
Q(  J-l )  sO. 

*ocj-d  =  i.oo7 
TCJ-1)  s  310.15 
VC  J-l )  S  o. 

XCJ)  s  X(J>1)«0X(J-1) 

10  CONTINUE 

Kl)  s  5.33E10 
PCI)  s  1 . 1 RL  ** 

P»IKT  *900 

WRITE  (6*510)  (J*  XCJ).  VCJ».  P(J)*  ROCJ)*  QCJ)*  T€J»*  O'lJli 
•EC J) •  J  =  l*  AO) 

CYCLE  si. 

00  no  CYCLE  si*  2000 
calculation  of  time  step 

OELT  s  50 C • 

D*  20  J  =  1*  AO 

ANOHCJ)  s  PC  J1MPC  J)*Q(  J)>‘GAHPRM-(RO(J>**?*2.3Af  10*CD.5ll« 

•  CROC  J)-l.  007)  -1.00  7)/CROCJ) -0.511  *CR0C J)-1.007))»«3) 

SPEEOCJ)  =  SORTCANUMC J)/ROCJ)) 

OTCJ)  =  FA ACT •CX(J*l)-XCJ) )/ SPEEOCJ) 

20  IF  COT C J)  «LT.  OELT)  OELT  =  OTCJ) 

CALCULATION  UF  CELL  BOUNDARY  ACCELERATION  AND  NEW  VELOCITY  VALUES 
WHILE  RETAINING  OLD  POSITION  AND  VELOCITY  VALUES 

UOOTCI)  s  0. 

DO  30  J  s  1*  AO 

30  UQOTCJ+1)  =  -A*CPCJ*mQCJM)-PC  J)-QCJ))/(CXCJ*2)-XCJ)>* 

•CROC J*l)*ROCJ)>) 

DO  AO  J  =  I*  A5 

OLDVCJ).  s  WCJ) 

AO  OLOXCJ)  s  XCJ) 

00  50  J  s  1*  AO 

50  VCJ*1)  s  VC  J^mUDOTC  J»l  )*OELT 

CALCULATION  OF  NEW  CELL  BOUNDARIES*  COMPRESSION  FACTOR  F*  OENSITY* 
ARTIFICIAL  VISCOSITY*  AND  PRESSURE  VALUES 

0.3  GO  J  s  1*  AO 

XCJ*1)  =  XCJ*l)t(COlOVCj*l)»VCJ»l))/2)*OCLT 
F  s  CCCL0KCJ*1)**AL»HA-C3LDXCJ))«*ALPHA)/CCKCJ*1))**ALPHA 
••CXCJ))** ALPHA)) 


VL- V»j 


—S.  **)&.*•**%>*  ^..V^V.^**'*-.  -  -.  *.•*.,**  ■•**  .V  fV^CfV'S.^ 


FP-M  a  F 

IF  CFPRH  .LE.  1  •  )  F**N  a  1* 

R3CJI  =  PO  <  J  >*F 

Q(J>  =  A**2*^0(J)»((X(JM)*X(JM**2l*(0CLT**(-2){ 

**lll.-fpm»»»2i 

iC  PfJI  =  PC J>-ANUM<J)*Cl.-FJ 


CALCULATION  OF  NEK  ENERGY  AMD  TEMPERATURE  VALUES 


0  70  J  a  1.  40 

£ H (  J)  =  C2. 30E10-2.32I10/OC  J) >« CO .4956*0. 5il/°0<J» >»•<-?> 
EMC  J)  =  C.5*CCl.«>lE5-l.32E5/ROCJII/C0.4e56*0.5U/ROCJIII**2 
EX  J)  =  C4PC J)-PHCJ))/0.51i*POCJ)l*EHCJ) 

TCJ)  a  T  ( J I  ♦  EC J)/1.593ie 
TIME  a  TIM£*0ELT 


PRINT  ROUTINE  FOR  RESULT S 

IF  CMGDCCVCLE •  1C)  .ME.  9.)  80  TO  100 
MRITE  C6*a00>  CYCLE*  TINE 

0:  FORMAT  C 1 X *  •PUN  =  ■*  T7*  14*  4X«  "TIME  =  •*  IPE9.3) 

PRIN“  $03 

OC  FOPHAT  CT10*  »X«*  T20*  «V«*  T30*  -P**  T40*  *PO«*  TSO*  «0«*  T60 
•  •T«,  T70*  "DT*  •  TOO*  •£•) 

W® ITE  C6#910  >  (J*  XCJ)*  VCJ>*  »CJ)»  RSCJI*  OCJ>*  TCJ)*  DTCJ)« 
•ECJ).  J  a  1*  40 
10  FORMAL  C1X*  13*  1P8E10.2) 

IOC  continue 
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Figure  E-l.  Sample  Output  from  EYESHOC 
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APPENDIX  F 


In  order  to  determine  the  validity  of  the  code  EYESHOC, 
a  point  explosion  was  modeled  using  the  code  and  was  compared 
to  the  results  obtained  by  Sedov  (Ref  25:238-248). 

Sedov  defines  his  shock  relations  in  terms  of  ratios  of 
pressure ,  density ,  radius ,  or  velocity  parameters ,  such  that 
the  comparison  of  all  shock  relations  can  be  made  using  these 
dimensionless  variables.  In  the  validation  of  EYESHOC,  the 
ratios  used  were  p2/p^,  $2^1*  an<*  r2^r°*  These  are  defined 
according  to 

P2  is  the  maximum  pressure  of  the  shock 
is  the  ambient  pressure  of  the  medium 
P2  is  the  maximum  density  of  the  shock 
is  the  ambient  density  of  the  medium 
r2  is  the  radial  distance  to  the  point  of  maximum 
pressure  of  the  shock 

and 

r°  -  <E0/Pl)1'3  (F-l) 

where 


Eq  is  the  energy  of  the  pulse  in  ergs 


To  validate  EYESHOC  based  on  the  dimensionless  variables 
introduced  herein,  a  point  explosion  was  modeled  for  an  ideal 
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gas  with  y  “  1.4.  This  was  used  rather  than  the  Griineisen 
equation  of  state  and  y  -  1.022  so  that  the  results  of  the 
comparison  would  be  apparent  and  not  hidden  because  of  the 
differences  in  the  conditions  of  the  test.  Figures  F-l  and 
F-2  present  the  results  of  this  comparison.  The  comparison 
was  made  for  an  energy  Eq  =  7.19  x  102°  erg  and  a  value  of 
r°  *  1.1  x  10s  cm.  As  is  noted  in  both  cases,  examination 
of  the  pressure  drop  across  the  front  and  the  density  behind 
the  front  as  a  function  of  shock  wave  radius,  the  dashed  line 
and  all  points  plotted  are  EYESHOC  results,  while  Sedov 
predictions  are  indicated  by  the  solid  line.  For  much  of 
the  range  covered  in  both  cases,  it  is  virtually  impossible 
to  distinguish  between  the  two  lines,  indicating  excellent 
agreement  between  EYESHOC  results  and  the  theory  presented 
by  Sedov. 

Based  on  this  comparison,  it  is  obvious  that  EYESHOC 
is  valid  in  its  treatment  of  a  point  explosion.  If  the 
assumptions  of  Chapter  IV  are  valid,  then  the  code  can  be 
assumed  to  be  usable  for  its  intended  purpose. 
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Radius  for  EYESHOC  and  Sedov 
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